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The heavens declare the glory of God; 
and the firmament sheweth his handywork. 
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their voice is not heard. 

Their line is gone out through all the earth, 
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Introduction 



Primordial Big Bang Nucleosynthesis (BBN) represents one of the greatest successes 

of the hot Big Bang model, along with the Hubble expansion and the Cosmic Microwave 
Background Radiation. Of the three, BBN probes the Universe to the earliest times, from 
a fraction of a second to thousands of seconds from its born, and its formulation predicted 
the existence of cosmic microwave background radiation Its emergence as a cosmo- 
logical cornerstone relies on the basic consistency of the predictions for the abundances 
of the light elements, such as L>, ^He, ^He, Li, with their measured abundances, which 
span over more than nine orders of magnitude. 

BBN took place in the early Universe when the temperature scale was less than 1 MeV. 
The key events leading to the synthesis of the light nuclides followed from the period when 
the weak interaction rates were in equilibrium, thus fixing the ratio of number densities 
of neutrons to protons; at temperatures T ^ 1 MeV this ratio was n/p ~ 1. As the tem- 
perature fell and approached the point where the weak interaction rates were no longer 
fast enough to maintain equilibrium, the neutron to proton ratio was given approximately 
by the Boltzmann factor, {n/p) ~ exp{— Am/T}, where Am is the neutron-proton mass 
difference. 

The nucleosynthesis chain begins with the formation of deuterium through the process 
p + n D + J. However, because of the large number of photons relative to nucleons 
(of the order of 10^''), deuterium production is delayed past until the temperature falls 
well below the deuterium binding energy of 2.2 MeV (the average photon energy in a 
blackbody is < >~ 2.7 T). 

The dominant product of BBN is '^He, resulting in a mass abundance close to 25 %. 
Smaller amounts of other light elements are produced: D and ^He at the level of about 
10^^, and ^Li at the level of 10^^*^ per hydrogen nucleus. In the standard model (see 
below) the abundances depend only on one free parameter, the baryon to photon ratio 
7], and remarkablya single value for rj ~ lO"^'^ may accommodate all observed orders of 
magnitude for ^He, D, ^He and ^Li data. Furthermore, the primeval yield of "^He is 
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also relatively insensitive to this quantity (pinning down r/ to 20 % pegs its value to 1 % 
precision) . 

The resulting abundances of the light elements are shown in Figure [l|, where we report 
the predictions for the '^He mass fraction Y, and D, ^He, Li abundances relative to H 
as function of rj. We also report the present observational situation (a complete discus- 
sion is presented in chapter ^). The general agreement between experimental data and 
expectations emerging from this figure, ranging over many orders of magnitude, is the 
great success of the Big Bang model which we referred to above. The essential problem 
in attempting to compare [0, [l^ the theoretical predictions with the observational data 
is that the primordial abundances have been significantly altered during the lifetime of 
the Universe through nuclear processing in stars and other galactic chemical evolution 
effects (we will discuss this subject in chapter §). The most stable nucleus, ^-fTe, grows in 
abundance with time since it is always created in stars, while D, the most weakly bound, 
is always destroyed. The history of ^He and ^Li is more complicated since these elements 
may be both destroyed and created during stellar evolutiom. To avoid corrections which 
are difficult to treat quantitatively, it is necessary to measure abundances in the most pri- 
mordial available material. Observations of light element abundances have dramatically 



improved over the past few years [10|. Although D and ^He abundances have about a 10 
% uncertainty and ^ Li data are even more uncertain, as an example we mention the fact 



that He data are now reaching a precision |11] of one per mille. 



This perspective, however, the fact that BBN is now entering in his maturity and precision 
era, demands similar improvements in the precision of the theoretical analysis, in order 
to reduce as much as possible all uncertainties in the predictions. An increasing precision 
in the measurements of '^He mass fraction at the level of 10~^ requires, for example, a 
reliability of neutron-proton conversion rates at the same level of precision of all other 
effects which are relevant for the neutron to proton ratio at the onset of nucleosynthesis 

This is the main goal of this thesis. 

The physics of BBN is well understood |T^: basically, it can be seen as a nuclear 
reactor in an expanding box. Therefore nuclear and particle physics and cosmology are 
the two basic inputs for studying primordial nucleosynthesis. Since BBN involves events 



0.26 




1 10 

Figure 1: Predicted abundances (solid curves) of primordial ^He, D, ^He and '^Li along 
with their experimental determination (rectangles) as a function of 7710 = lO^^rj. The two 
curves for '^He correspond to the la experimental values for the neutron lifetime, while 
the theoretical uncertainties on "^Li are as discussed in Ref. Uncertainties on D and 
^He curves cannot be appreciated on the scale of this figure. The filled zone is the range of 
T] representing the agreement between BBN theory and measurements (taken from fl^ ). 
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that occurred at temperatures of order 1 MeV, it naturally plays a key role in forging the 
connection between cosmology and nuclear and particle physics. It is interesting to note 
how the increasing interaction between particle physics and cosmology has largely resulted 
from the establishments of "standard models" in both fields which satisfactorily describe 
all known phenomena. Regarding cosmology, there are of course many possibilities for de- 
partures from its standard model Q, e.g. an inhomogeneous nucleon distribution or non 
zero neutrino chemical potentials. However, recent developments do not motivate such non 
standard scenarios and moreover they are highly constrained by the observational data. 
It is therefore reasonable, and we hereafter shall do so, to adopt the standard picture. In 
our calculations we will also use the standard model of strong and electroweak interactions 
[p!^, but it is very intriguing to observe how BBN can constrain new physics beyond the 
standard SU{3)x SU{2)xU{l) model 0. This is well illustrated by the BBN limit to 



the number of light neutrino species |17]. The amount of synthesized "^He strongly de- 
pends upon the expansion rate: a faster expansion rate leads to the production of a larger 
amount of '^He since nuclear reactions begin earlier when the neutron fraction is higher. 
The expansion rate itself is determined by the energy density of relativistic particles and 
the larger the latter, the larger the former. Now, since the energy density increases by 
adding new neutrino species, this can overproduce "^He, thus violating the observational 
evidences. The BBN limit on the number N,y of light neutrino flavours is close to three, 
although more work, both observational and theoretical, is still needed. It is impressive 
and interesting to compare this result with the LEP limit based upon the shape of the 
Z° resonance: = 2.994±0.012 @. While it is unhke that the BBN limit will ever 
achieve such precision, the cosmological and laboratory limits are complementary. The 
neutrino limit based upon the shape of the counts the number of "active" particle 
species that have a mass smaller than half the mass, weighted by their coupling to the 
Z^. Differently, BBN constrains the energy density due to any relativistic particle specie 
around the time of primordial nucleosynthesis and thus is sensitive to any particle species 
lighter than about lAIeV. 

Other several interesting constraints on physics beyond the standard model of elementary 
particles, coming from BBN considerations, can be obtained but this is not our subject. 
We only stress the fact that accurate BBN predictions, along with precise experimental 
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observations, greatly helps also our understanding of fundamental physics. 

Summarizing, this thesis is devoted to the study of precision effects intervening during 
primordial nucleosynthesis. Our present goal is to give an estimate of the primeval ^He 
mass fraction confident, at least, at the third significant digit (the actual experimental 



accuracy). The following step is to write down a new numerical BBN code [19| employing 
our results to further improve light element abundance predictions. 

The thesis is organized as follows. Chapters 1 and 2 review standard cosmology and, 
in particular, primordial nucleosynthesis. The subsequent three chapters are devoted to 
the calculations of the mentioned precision effects. In particular we discuss the calculation 
of QED radiative corrected Born rates for the weak reactions fixing the neutron to proton 
ratio at freeze out, along with the calculation of corrections to this quantities arising from 
finite nucleon mass and QED thermal plasma effects. Finally, in chapter 6 we summarize 
our results and give our prediction for the '^He mass fraction. 

Throughout this thesis, unless otherwise specificated, we use natural units in which h = 
c = k = l. 
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Chapter 1 

The standard cosmology: an 
overview 



The standard cosmological model is based on three observational pillars: 

i) the uniform distribution of matter in the Universe on large scales and the isotropic 
expansion of it that maintains the uniformity; 

ii) the existence of a nearly uniform and accurately thermal cosmic background radia- 
tion (CBR); 

iii) the abundances (relative to hydrogen) of the light elements D, ^He, ^He, Li. 

As we will see, the validity of the Bubble's expansion law, namely the proportionality 
between the observed red-shift z (see below) and the distance d from the Earth of a given 
source Q 

H^d z . (1.1) 

through the Hubble constant Hq, out to red-shifts z ~ 0.2, supports the general notion 
of an expanding Universe [^0|. The Hubble constant Hq, defining the present expansion 
rate, is usually parametrized as 

Hq = imho Kms-^ Mpc-^ , (1.2) 

where, observationally |21], |22| 

0.6 < ho < 0.8 . (1.3) 

^In astronomy, this distance is defined as d? = £/(47rjr), with C the absolute luminosity (energy per 
unit time produced by the source in its rest frame) and T the measured energy flux (energy per unit time 
and surface measured by a detector in the expanding Universe). 



1 



The distribution of matter and radiation in such a Universe is observed to be homogeneous 
and isotropic when averaged on scales exceeding a few hundred Mpc Instead, the 

CBR provides a firm evidence of a hot, dense beginning of the Universe itself ("Big Bang"); 
the spectrum of the CBR is a perfect blackbody, at temperature [E4| 



To = 2.728 ± 0.002 K , (1.4) 
with deviations that are less than 0.03 %. 

Observations of the primordial light element abundances finally provide the complete 
success of the hot Big Bang model, and they will be discussed in the next chapter. 

1.1 Homogeneity and Isotropy 

Space-time events in a homogeneous and isotropic Universe are described by the maximally 
symmetric Robertson- Walker metric which, in the comoving reference frame (defined by 
the property that an observer at rest in this frame has constant spatial coordinates in 



time), takes the form [25] 



ds'^ = - R^{t) L _ + r'^dO'^ + sin^^^M , (1.5) 

where t G [0, -|-oo) is the proper time measured by an observer at rest in the comoving 
frame, and with r G [0, 1], £ [0, vr], E [0, 27r] the (dimensionless) spherical coordinates 
in this frame. The parameter k = -|-1,0,-1 gives the spatial curvature of the Universe. The 
homogeneity and isotropy of space allows us to describe the dynamics in the space-time 
in terms of only a quantity (with dimension of a length), the cosmic scale factor R{t), 
depending on the time t only. 

Since, in the Robertson- Walker metric, R{t) is a function of time, the distance between 
two space points depends on time; if the space described by this metric is expanding or 
contracting, the motion of a particle is influenced by this expansion/contraction. In fact, 
let us consider, for example, a particle in free motion with momentum pi at time ti; 
since the momentum has dimension of an inverse of length, at a later time t2 the particle 
momentum will be rescaled by a factor proportional to R^^, that is 
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Furthermore, due to the finite propagation velocity of hght signals, in an expanding/ contracting 
Universe it is meaningful to introduce the concept of a "distance to the horizon" dn (t) : for 
a comoving observer, this is the distance at which a light signal emitted at t = reaches 
him at (or before) time t. If dnit) is finite, then there are sources from which light has not 
yet reached us [^]: a boundary exists (termed "horizon") between the visible Universe 
and the part of Universe from which light signals have not reached us. The distance dnit) 
can be calculated in the following way: from the isotropy of space, we can consider an 
observer at r = for which d9 = dcj) = 0; hence 

dH{t) = Va^dr = R{t) _ . (1.7) 

For a light signal ds^ = 0, thus 

and correspondingly 

The explicit form of dnit) depends on the expression for R{t) which is determined by the 
cosmological equations (see next section); however, we note that the behaviour of R{t) 
near the initial singularity determines the finiteness of dnit). In the standard cosmology 
this is indeed finite and dnit) ~ t. 

Another interesting property of an expanding/contracting Universe is that if at a given 
point Pi a photon is emitted with a wavelength Ai, at a distant point Pq it will be detected 
with a different wavelength Aq. In a space described by the Robertson- Walker metric, on 
dimensional grounds we have in fact that 

^ - fi ■ 

The red-shift (or blue-shift) z of a given object is defined as the ratio between the variation 
in the detected wavelength and the emitted wavelength: 

Ao — Ai 



Ai 



(1.11) 



From (1.10) we then have 



An increase (or decrease) in the cosmic scale factor then leads to a red-shift (or blue-shift) 
of the light from distant sources. The Hubble law ( |1.1[ ), relating the distance to a galaxy 
with the observed red-shift through the Hubble parameter at the present epoch 

- fi • 

to being the present time, is a direct consequence of the Robertson- Walker metric (its 
derivation is reported in for example); it is only an approximate relation ||2^, and 
corrections are necessary for cosmological large distances. The present age is set by the 
Hubble time 

~ 9.778x10^ /iQ^r , (1.14) 
corresponding to a local spatial scale for the Universe 

Hq^ ~ 3000 Hq^ Mpc (1.15) 

(Hubble radius). 

1.2 Cosmological equations 

As anticipated in the previous section, the cosmic scale factor governing the dynamics of 
the processes in the Universe is determined by the Einstein field equations relating the 
energy-momentum tensor T^i, to the space-time curvature: 

Rfiu - ^RcQ/iu = ^irGT^.y . (1.16) 

Here G = Mp'^ is the Newton gravitational constant, while -R^i, is the Ricci tensor and 
Rc its trace (the scalar curvature). For a perfect fluid, as the Universe is assumed to be 
(see the next section), the energy-momentum tensor takes the form 

Tf,u = -Pgfiu + (P + P) UfiUu , (1.17) 

p and p being the pressure and energy density respectively, while = dx^/ds is the fluid 
4-velocity. In a comoving frame (u^ = (1,0)), Eq. ( |1.17 ) simplifies to 



T 



( p \ 

p 

p 

V p / 



(1.18) 



Note that from homogeneity and isotropy of space it follows that p and p, as well as R, 
only depend on the time t. The Ricci tensor for the Robertson- Walker metric is instead 
given by 



Rij 



R 







' R R^ k ^ 

R ^ '^R^ ^ ^R^ 



9ij 



and 



Rr 



k 



jR R^ 
\ R ^ R2 ^ R^ 



(1.19) 
(1.20) 

(1.21) 



Substituting these expressions in the Einstein equations ( 1.16| ), we then obtain the follow- 
ing two independent equations 0: 



R' 
R 



+ 



R 
R 



i?2 
AttG 



SttG 
{p + 3p) 



(1.23) 



(1.24) 



For both matter and radiation p + 3p is positive, thus R is always negative. Since we 
know that now i? > (the Universe is now expanding) , at a given remote (but finite) time 
t^B in the past we have had that RitBs) = 0; this event is called the "Big Bang", and it 
is usually chosen as the starting reference time {tsB = 0). 

In the Eqs. ( |1. 23 ), (1.24) the energy-momentum conservation equation DyT^'^ = {Dy 
being the covariant derivative in the Robertson- Walker metric) is also contained; in fact, 
deriving ( |1.23| ) with respect to time and substituting ( |1.24D , we get 



p = -3H (p + p) , 
where we have introduced the Hubble parameter (depending on time) 

H = S^ . 

R 



An alternative form of (1.25) is as follows: 



d(^pR 



P 



d{R^) , 



The equation ( 1.23 ) can also be cast in tlie intriguing form 

GM _ fc 
" 2 



1 2 



R 



(1.25) 



(1.26) 



:i.27) 



(1.22) 



with V — R and M — ^-kR? p, showing that the sign of k is deciding about the final destiny of our Universe. 



Radiation 


P = 




P ^ 




Matter 


P = 





P ^ 




Vacuum 


P = 


-P 


P ^ 


-' constant 



Table 1.1: Equations of state for radiation, matter and vacuum. 

expressing just the first law of thermodynamics for the expanding Universe. 
Equations (|L2|) and (|L25|) 

SvrG k 



P 



3 ^ 

3H {p + p) 



(1.28) 
(1.29) 



are usually considered as the two basic independent equations governing the dynamics 
of the Universe. They are known as the cosmological Friedmann-Lemaitre equations. 
Obviously, only two equations are not sufficient to determine the three unknowns R{t), 
p{t), and p{t). The third relation to add to the Friedmann-Lemaitre equations is the 
equation of state relating pressure and energy density, that can be written as 



p(t) = wp{t) , 



(1.30) 



where, for simplicity, we will assume the coefficient w to be time independent. From ( 1.25| ) 
we then have 

p ~ i?-3 (!+«') . (1.31) 

The interesting case of radiation (relativistic particles), matter (non relativistic particles) 
and vacuum are summarized in Table |L1] . 

1.2.1 Evolution of the Universe 



The evolution of the Universe is determined by the curvature term k/R"^ in ( |1.28 ) which 
is positive, zero or negative if the energy density is greater than, equal to or less than the 
critical density 

Pc = ^ , (1.32) 



SttG 
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respectively, since Eq. (|1.28| ) can equivalently be written as 

k 



with 



O - 1 , (1.33) 



Vt = . (1.34) 

Pc 

For k = —1, E? is always strictly positive and R ^ t as t ^ oo. Instead, for /c = 0, 
B? goes to zero as i? — > oo while, for A; = — 1, i?^ drops to zero at Rmax = \/3/ (SirGp) 
after which R begins decreasing. Thus Q < 1 corresponds to an open Universe which will 
expand forever, = 1 is a flat Universe which will asymptotically expand to infinity while 
> 1 corresponds to a closed Universe which will eventually recollapse. 
The critical density today is (from ( |1.32| )) 

, 4 



PcO 



(2.999x10^^2 V/io Gey) = 1.054x 10"^ /ig GeF cm"^ . (1.35) 



From dynamical measurements of the present energy density in all gravitating matter (and 
excluding the nowadays negligible contribution of relativistic particles) we deduce |2^, ^ 

0.1 ^ 1 . (1.36) 

The present energy of photon background alone is known with a very good accuracy from 
measurements of the temperature of the cosmic microwave background radiation (|1.4|), 



which gives (see sect. 1.3.1 ) 



and then 



,,„ = ^;|l,2.02xlO-(J^)'G.y' , (1.37) 



= ^ . 249X10-. . (1.38) 

Including the contribution of a primordial background of three massless neutrinos (see 



sect. 1.4.1) we have 



= n^o + ^.0 - i.Gsn^o ~ 4.18x10-5 (y^^ ■ ^^-^^^ 



Comparing (ll.36| ) with ( 1.39 ) we then see that the present Universe is dominated by non 



relativistic particles; however, as we will see in the next subsection, the evolution to this 
situation has been highly not trivial. 
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1.2.2 Driving the expansion 



The equation (1.28) for the expansion rate of the Universe can be conveniently normalized 
to present values of the quantities involved {Rq = R{to), Hq = H{to) and = p/pco, to 
being the present time), obtaining 

^ "o(., - ^ (f . (1.40, 



From Table 1.1 we then see that 



Qoit) = nnoit) [^j + ^Moit) (^-^j (1.41) 

(the pedices i = R,M indicating radiation and matter respectively), where themselves 
are in general functions of time which can be assumed, for simplicity, to be step function 
-hke: 

Qioit) = I . (1.42) 

Thus we get 

(^)^ - (f )* . (f )' - 4, i^r . 

Given the different powers of R for the terms in ( |1.43| ) determining the expansion rate, we 
see that the early Universe (R « Rq) was dominated by radiation (RD era), afterwards 
matter becomes dominating (MD era). The transition between the RD and MD era can 
be approximately dated back at Req ~ 10~^i?o- Since in the first stages of evolution of 
the Universe the curvature term can be neglected, we find that 

RDera ~ (1.44) 

MDera ~ R'^ (1.45) 



and thus, from (|1.33|) , 



= {l + z)~^ MDera 



|0 - 1| ~ <( . (1.46) 

^(A)' = 10^(1 + -)-^ RDera 
Note that, apparently, at very earlier epochs the Universe was very nearly critical (of the 
order of one part over 10^). However, at these times a period of "inflation" (exponential 
growth of the expansion) is believed to have occurred [^9| , but a thorough analysis of this 
issue is beyond the scope of this thesis. 
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1.3 Thermodynamics of the Universe 



In this section we will study the properties of the Universe considered as a thermody- 
namic system composed by different species (electrons, photons, neutrinos, nucleons, etc.) 
which, in the early phases, were to a good approximation in thermodynamic equilibrium, 
established through rapid interactions. Obviously, coming back to the past, decreasing the 
cosmic scale factor we have an increase of the temperature. In this thesis we will mainly 
concern ourselves with cosmological processes occurred during the RD era, thus the dis- 
cussion is greatly simplified by assuming the Universe as an ideal gas. This assumption is 
justified by the fact that, as we will see in the next chapter, the particle densities do not 
usually become high enough for many-body interactions to be important. 

1.3.1 Equilibrium thermodynamics 

In a gas of a given specie with gi internal (spin) degrees of freedom and energy Ei = 



+ m?, kinetic equilibrium is established by sufficiently rapid elastic scattering pro- 
cesses; in this case, for an ideal gas, the equilibrium phase-space density is 

hip) = M^^^] ilV' , (1.47) 



Ti 

where +/- refers to Fermi-Dirac/Bose-Einstein statistics and is the chemical potential. 
In general each specie has its own equilibrium temperature Tj, and the entire Universe 
can be represented as a plasma with different temperatures. However, if several species 
strongly interact among them, they will reach a mutual chemical equilibrium and a com- 
mon temperature; this is indeed the situation at early times. As the Universe expands 
and cools down, some species may start interacting more and more weakly and eventually 



decouple. As we will see in sect. 1.4, we can consider the photon temperature as the 
plasma reference temperature T of the Universe. 

In chemical equilibrium, established by processes which can create and destroy particles 
(differently from kinetic equilibrium), the chemical potential is additively conserved. So 
it is zero for particles such as photons and which can be emitted and absorbed in any 
number and consequently opposite for a particle and its antiparticle which can annihilate 
into such bosons. 

The quantities of interest are the number density, energy density and pressure of a given 
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specie, defined in general as: 



= gi I -^Mp) , (1-48) 



9: 



(27r)= 

cPp IpI^ 



fiip) ■ (1-50) 



(27r)3 3Ei 

In kinetic equilibrium the phase-space density fi [p) is given in ( |1.47D , and these quantities 
evolve according to temperature. 

For non relativistic species (Tj ^ rrii) we have (for both Fermi-Dirac and Bose-Einstein 
statistics) 

3 

^ 5i [-1^ j e T , (1.51) 
Pi ~ m (nii + It,] , (1.52) 



2 

Pi ~ riiTi <^ Pi (1.53) 

(and so we recover Boltzmann statistics). The average energy per particle < Ei >= pi/rii 
and net number density are instead given by 

<Ei> ~ + ^Ti , (1.54) 
/ miTA i . , Hi _iiH 

(assuming = = /ij) . 

For relativistic species (Tj ^ rrij) we obtain: 



non degenerate case (Tj /ij 



3 1 
4 



nj ~ (1.56) 

I m9^Tt ED 

p., { (1.57) 
m9^Tf BE 
1 

Pi ^ 3 A (1-58) 

r I f^C(3)Ti ~ 3.15 Ti FZ? 

< Fi > ~ <^ (1.59) 

[ MC(3)ri ~ 2.70 Ti BE 

nj+ — Tii- ~ (1.60) 
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weak degeneracy (/ij < 0, < Tj) ^ 



n^ ~ \giTfe^^^ (1.61) 

p, ~ 2,\giTte^^ (1.62) 

Pi ^Pi (1-63) 

<^i> ~ STi (1.64) 

ni+ - m- ~ 2 ^ fifj sinh ^ (1.65) 



degenerate case (Tj <C ;Uj 



Pi - -^9if4 (1-67) 
1 

Pi ^ 2 Pi (1-68) 

<Ei> ~ ^;Ui (1.69) 
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In the relations above C,{x) is the Riemann zeta-function, and ^(3) ~ 1.2021.... 

Another important quantity for the evolution of a thermodynamic system is its entropy 
S which, for thermal equilibrium, is defined by the second law of thermodynamics: 

TdS = 6Q . (1.71) 

Assuming zero chemical potentials, for a given comoving volume element {V = R^) we 
have 

TidSi = d{piV) + VidV = d{{pi + Vi)V) - Vdvi (1.72) 

(note that pi and pi are the energy density and pressure of species at equilibrium). The 
entropy function is subject to the integrability condition 



dTi dV dV dTi ' 



(1.73) 



■^For bosons with > Q Bose condensation may take place 
''The foUowing resuhs hold for both FD and BE 
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which exphcitly imphes that 

Pi + Pi ^ 1 dpi 1 dpi 1 dpi 



or 



dp.. 



Ti dTi Ti dTi Ti dTi 
Pi + Pi 



dTi . 



(1.74) 



(1.75) 



This relation is in fact verified by ( 1.49] ), ( |1.5[1| ) with the equihbrium distribution ( |1.47| ). 
Substituting in ( 1.72|) we then have 

dS. = d(v^-^P^^ 



(1.76) 

from which we found that the entropy per comoving volume {V oc B?) is defined, up to 
an additive constant, by 



Si = V 



Pi + Pi 



R 



Pi + Pi 



(1.77) 



(1.78) 



Note that from the first law of thermodynamics ( 1.27 ) and (|1.72| ) 

dSi = , 

i.e. in thermal equilibrium ^ the entropy in the comoving volume is conserved. 

It is also useful to define the entropy density in the comoving volume, 

^ 5i _ Pi + Pi 
- V ~ Ti 

which, from ( |1.78| ), varies as 

S^ OC i?"^ . (1.80) 

The total number of particles of a given specie in the comoving volume Ni = Vui = E?ni 
is thus proportional to rii/si, 



(1.79) 



Ni oc 



(1.81) 



Finally, if chemical potentials are not zero, Eq. ( 1.71[ ) specializes into 

TidSi = d{piV) + pidV - pid{niV) , (1.82) 
and then the entropy is now given by 

J-i 

which replaces Eq. (1.77). 

^During phase transitions entropy is not conserved, in general. 
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1.3.2 Energy density and pressure in the Universe 

In the Friedmann-Lemaitre equations ( |1.28| ), ( |1.29| ), assuming a perfect fluid, the total 
energy density p and pressure p of the cosmological fluid are involved; let us evaluate 
these quantities considering all the particles in thermal equilibrium in the Universe and 
express them in terms of photon temperature = T: 

where the sums run over all species and Ui = Ei/Ti, Xi = rrii/Ti, yi = pi/Ti. Note that 
from ( 1.51| ), ( |1.52 ), non relativistic particles contribute negligibly to the energy density in 



the radiation dominated era, since their energy density is exponentially suppressed with 
respect to the case of relativistic particles; thus we can neglect the contribution of non 
relativistic species in the sums above. Assuming all species non degenerate, we then get 



, (1.86) 



where in contribute only the relativistic degrees of freedom 

[B = bosons, F = fermions). Note that g^, is in general a function of T, since the number 
of degrees of freedom becoming relativistic at a given temperature depends on T itself. 
Moreover, at a given time, not all (relativistic) particles in the bath are in equilibrium at 
a common temperature T. A particle will be in kinetic equilibrium with the background 
plasma (that is Tj = T) only as long as its interaction with the plasma is fast enough; 
although the conditions for this to occur will be discussed in sect. |1.4[ , it is obvious that 
these involve a comparison between the particle interaction and the expansion rate H. 

1.3.3 Total entropy density in the comoving volume 



From (IlTQI) and ([Ts^ ), (|L87|) we can evaluate the total entropy density in the comoving 



volume, obtaining 



27r2 
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with 

9*s - E5^(§)' + ^E5^(§)' ' (1-90) 

and the sums again run only over the relativistic degrees of freedom in equihbrium (in 
the considered approximation). Note that s in ( 1.8S| ) can be parametrized in terms of the 
photon number density as follow 

s = 45(^(3-) 9*s n^i ^ l.^Qg^sn^ ■ (1-91) 

From entropy conservation we can now obtain the scaling law relating the cosmic scale 
factor with the temperature; in fact from ( |1.78| ) and ( |1.89D we get 

T ~ g^s^'^R-^ . (1.92) 

We stress that if g^g depends on T, then the adiabatic invariant RT = const, for pure 
expansion does not hold. 

1.3.4 Decoupling 

Non equilibrium phenomena are properly described by the Boltzmann equation for the 
phase space distribution functions of the given species involved in them. Such description 
is reported in appendix Here we make a very simple and general analysis of non 
equilibrium phenomena occurring during the evolution of the Universe. 
A specie i is said to decouple at a temperature (corresponding to a given R^ and time 
to) when it goes out of equilibrium with the background plasma. The conditions for this 



to occur will be considered in sect. 1.4. Of course no particle is ever truly decoupled 
since there are always some residual interactions with the bath, but it is a convenient 
assumption to consider that a specie decouples at a given temperature T|). 
Let us first consider species that are relativistic at decoupling (mj < T}-,) Q In this case 
the distribution at decoupling results to be the equilibrium one: 

Subsequently the decoupled i particles will expand freely without interactions, hence their 
number in a comoving volume is conserved and their thermodynamical quantities are 



''For simplicity we consider only the non degenerate case 
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functions of the scale factor R alone. Although non interacting, after decoupling their 
phase space distribution will retain the equilibrium form as long as the particles remain 
relativistic, since from (|1.6D we have 



mo) = Ei (^A^ , (1.94) 



and thus for t > 



/.(p,T.<Ti)= (e.p{^j±l) =(„p||}±l)-' , (1.95) 

This is because of the adiabatic invariant (for pure expansion) 

T^ = Th (^) ~ R-' . (1.96) 

Initially the temperature Tj exactly follows the photon temperature but, as the Universe 
cools below some mass thresholds, the corresponding massive particles will become non 
relativistic and annihilate. This will heat the photons and other interacting particles but 
not the decoupled i particles, so that Tj will drop below T and consequently rii/n^ will 
decrease below its value at decoupling. 

For species that are non relativistic at decoupling (m^ » T^), ( |1.93| ) continues to be valid 
and after decoupling the phase space distribution will again retain its equilibrium form, 
but now Ei^tf)) ~ p\,/2mi and from (|1.6|) we have 

E^{tD) = Ei (^-^y , (1.97) 



thus ioi t>tD 



EiR 1 \ / f E., 



where 



T^ = Th[^)' ■ (1-99) 



From the discussion above it is then clear that ultrarelativistic and non relativistic species 
retain their equilibrium distribution after decoupling with temperatures varying as in 
( |L96| ) or ( IQgl) . This is no longer true for species with intermediate energy at decoupling, 
Tjj ^ mi. In fact for these we have 



ElM. = J_Jpl+m^ = J_./fA)%2 +^2 ^ El , (1.100) 
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which "defines" a momentum dependent temperature, indicating thus a non equihbrium 
distribution. 

1.3.5 Time - temperature relationship 

The useful relationship between the time t and the background (photon) temperature 
T in the radiation dominated era can be obtained straightforwardly by integrating the 
Friedmann-Lemaitre equations by means of entropy conservation. In fact, 

/■^(*) 1 dR , , 

' = /„ Mil ' <'-l°" 

by using the Eq. ( 1.2^ ) (neglecting the curvature term in the RD era) and (1.86), giving 



H 



I 47r3 
45M|, 



and entropy conservation ( 1.78 ) with Eq. ( |1.89| ) 



dR dT dq^ 
3^ = 3— + ^ 



we find 



t 



U5Mj, „i/2 
9* 



1 + i ^ '^9*s \ dT 
3 g^s dT J T^ 



(1.102) 



(1.103) 



(1.104) 



V 47r2 JO 

During the periods when both g^, and g*^ are approximately constant (i.e. away from 
phase transition and mass thresholds where relativistic degrees of freedom change) the 



_i/2 flMeVV 



sec 



relation (1.104) simplifies to 

t ~ 2.42 5* ' [-^jT- ) 
or, quite approximately, t{sec) ~ T^'^{MeV). 

1.4 Thermal evolution of the Universe 



(1.105) 



The cosmological model of Friedmann-Lemaitre, as described in the previous sections, 
does not allow a Universe in thermal equilibrium since, in this case, the Universe itself 
would be stationary, which is not the case. However, sufficiently away from the Big Bang 
event, we can nevertheless approximate the evolution of the Universe as made of several 
subsequent phases of (different) thermal equilibrium with temperature T mainly varying 
as R^^ . Thermal equilibrium is realized if the reactions between the particles in the heat 
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bath take place very rapidly compared with the expansion rate, set by H. Thus denoting 
with T = TLi < av > the thermal average of the given scattering reaction rate {a and v 
being the cross section and the particle velocity, respectively), the (approximate) condition 
for having equilibrium is given by |^2[ 

r > H . (1.106) 

Of course the opposite condition T < H is not a sufficient one for the exit from the equi- 
librium since, as we have seen, for example, a relativistic low interacting specie maintain 
its equilibrium distribution but with T ~ R^^. Nevertheless the decoupling temperature 
Td, as used above, is defined by 

r(Tz)) = H{Td) . (1.107) 

The correct way of proceeding is to solve the Boltzmann transport equations for the given 
specie (see Appendix^), but here we are interested in a semi-quantitative discussion, and 
the criteria above fit well. 

The scattering rate depends on the particular interaction experienced by the particles, 
while (in the RD era) the expansion rate is approximately 

H — . (1.108) 

For interactions mediated by massless gauge bosons (for example photons or W^,Z^ before 
the electroweak phase transition) there is no particular mass parameter, so that for reac- 
tions of the type a + b — > c + d the cross section is approximately given by < cr >~ q'^T"^, 
a being a dimensionless coupling constant. Thus, for example, for relativistic particles in 
the radiation dominated era (and n ~ T'^) we have 

r ~ a^T . (1.109) 

Hence for this type of interactions the decoupling temperature is set by Td ~ a'^Mp and 
for T ^ a^Mp equilibrium is established, while for T ^ a^Mp the reaction rates take place 
slowly (they are frozen out). 

For interactions mediated by gauge bosons with mass Mx we now have < a >~ G\T'^ ~ 
a'^MZ'^T'^ and 

T ~ . (1.110) 
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In this case the decouphng temperature results T/j ~ {Mj^/a^MpY/^ and for temperature 
above this (but below Mx) equilibrium is maintained, while for T < To the reactions are 
frozen out. 

In practice, apart from the massless neutrinos (whose decoupling will be studied be- 
low), all particles in the Standard Model of elementary particles are strongly coupled to 
the thermal plasma (of 6^,7) while they are relativistic. This is why we can reliably 
describe the evolution only in terms of the 'plasma' (photon) temperature T and neutrino 
temperature T^,. Thus the thermal history of the Universe, as pictorially reported in table 
[1.2| , can be fairly reconstructed reliably back to the Fermi scale but, with some cautions, 
also nearly up to the GUT scale. This happens because of the dilute RD plasma approxi- 
mation, in which non relativistic particles of plasma are (forced) in equilibrium with the 
relativistic species with negligible abundances. Some uncertainties arise, however, when 
the ideal gas approximation breaks down, that is at phase transitions associated with 
symmetry breaking and at very high temperatures, where equilibrium is not achieved. 

1.4.1 Neutrino decoupling 

The above discussion can be interestingly illustrated by the example of the decoupling of 
(massless) neutrinos. They are maintained in equilibrium by reaction such as vV — > e'^e~ , 
ve — > ve and so on, whose averaged cross sections are of the order of < ex >~ G'pT'^, 
where Gp is the Fermi coupling constant. From ( |1.107] ), ( |1.108| ), ( |1.110| ) we then find the 
neutrino decoupling temperature 

~ IMeV (1-111) 



(a more careful estimate of the cross sections [33, 34] gives T^'^ '^^ ~ 3.5 MeV and T^*^ ~ 



2.3 MeV, the higher decoupling temperature of ffi^fT depending on the fact that they 
interact only through weak neutral current, while for there is also a charged current 
contribution) . Thus above ~ 1 MeV neutrinos are in equilibrium with the plasma of 
photons with Ty = T and (from ( 1.56D ) riy = (3/4)n^, afterwards they decouple from 



the plasma and their temperature Ty scales approximately as R~^. Subsequently, as the 
temperature drops below ~ 0.5 Mel/, e"^ and e~ annihilate, heating the photons but not 
the decoupled neutrinos. In these three phases we have respectively 

gl^) = g.{T>T-^) = 2 + ^(4 + 3x2) = | , (1.112) 
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gf^ = 5*(rne<T<Ti^) = 2 + ^4 = y , (1.113) 
gf^ = g.{T<me) = 2 . (1.114) 

The evolution of the neutrino temperature through the period of annihilation can be 
easily computed by entropy conservation. Note that from ( 1.78| ) the entropy of each specie 



in equilibrium in the comoving volume is conserved but, in the present case, also the neu- 
trino entropy after their decoupling is conserved since, after decoupling, the total neutrino 
number in the comoving volume does not change anymore (assuming stable neutrinos). 
Thus from the conservation of the total entropy S = Sj + Sy and neutrino entropy Sy, 



it follows that also the entropy Si associated to the species still interacting |15|, given 
approximately by (see ( |1.91 ), ( 1.89| )) 



Si g^Ny g.R^T^ , (1.115) 

is conserved, i.e. 



^(2) ^ ^(3) ^ ^-^ -^^g^ 



implying, from ( |1.113|) , ( |1.114| ), 



(i?r)(3) = (^^y (RT)^^^ . (1.118) 

However, (RT)^^'' = {RT^)^'^^ = const. = {RT^f'^^ (since ~ R-^), for which, after 
annihilation 



, . , • (1.119) 

Note that after annihilation there are no other relativistic species that can become non 
relativistic (altering the effective degrees of freedom) so that ( |1.119| ) is just the relation 
between the present values of 7 and v temperature, as well as ( |1.117] ) relates the present 
number of 7 to that preceding annihilation. Furthermore, we also have from ( 1.117] ) 



\nr^j 11 \n^j 11 4 11 ^ ' 

From this value of T till present days, neutrinos remain relativistic and therefore continue 
to retain their equilibrium distribution; hence the degrees of freedom characterizing the 
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present energy density and entropy (7 + 82^) are given by 



4 

gl = 2 + ^ (3x2) (^^y 3.36 (1.121) 
als = 2 + ^ (3x2) (1) = (1.122) 



from which we deduce, for example, {T ~ 2.73 i^) 



nO = ^^T^ ~ 422 (1.123) 
2vr2 

4 = -^alsT^ ^ 2970 (1.124) 

PR = -^alT^ - 8.09x10^34^^^-3 (;l.125) 

n\hl ~ 4.31xl0~^ . (1.126) 



The present quantities related to neutrinos are obtained from 



3_ 

iT 

1.96A" . (1.128) 



nl = (1.127) 



A more careful analysis of the -T relation is conducted in appendix ^ where we report 
the generalization of Eq. ( 1.11S| ). 
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Table 1.2: History of the Universe. 



Time 


Tomporat urc 




Parliclc coiitcnil 


today 


T ~ 2.73 ii: 




7 + 3 decoupled v 


t ~ 10^ '' s 


T ~ 10"^ eF 


Galaxy formation 




t ~ 10^^ s 


r~ leF 


Matter-Radiation decoupling 
Atom formation 


5? 


t ~ 10^^ s 


r~ lOeV 


Pmatter ~ Pradiation 
Structure formation 


?? 


i ~ 10'' s 


T - 10^ ey 


Planck spectrum established 






T~ me 




add 




T ~ 1 MeF 


Light element nucleosynthesis 






T ~ 2 ^ 3 MeF 




vs become interacting 




T ~ 




add p^ 




T ~ 




add vr^, vr" 


t - 1(J-" .s 


150^ lOO.Vc V 


Quark-hadrou I rausil ion 










7,3i/, e±,/i± 
u, u, d, d,s,s, gluons 




T ~ mc 




add c, c 








add 




T nib 




add b, h 




T ~ miv,z 




add T^±,Z" 




T ^ mt 




add t, t 


t ~ 10-^^s 


300 GeV 


Electroweak phase transition 


add iJ" 


t ~ 10-^4 s 


T ~ 10^'^ GeF 


GUT symmetry breaking 
Inflation 


add X, Y 


t - 10-« s 


r ~ lO^y GeV 


Quantum gravity 


? 
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Chapter 2 

Primordial Nucleosynthesis 



2.1 Preliminaries 

During about the first 20 minutes in the evolution of the Universe, conditions were 
favourable for the synthesis of (significant abundances of) light nuclides, such as D, "^H, 
^He, ^He, '^Li, ... . The onset of nuclear reactions which build up the light nuclei takes 
place slightly after the decoupling of the weak interactions (see table |1.2| ) which keep neu- 
tron and proton in chemical equilibrium. This is a necessary step since, as long as free 
nucleons are in equilibrium, the (equilibrium) abundances of all bound nuclei are quite 
negligible, due to the very high value of the entropy per nucleon {s/un ~ lO^'^). 
Introducing the subject, we will start with some definitions and comments. The number 
density ua of a given nuclear specie {Z, A) has been defined in ( |1.48| ) and, since the specie 
is non relativistic, it is given by the expression ( p.. 51 ): 



riA ^ gA [ -T, — e T . (2.1) 



2vr 

The baryon number (density) is the total number of nucleons (bound or free) in the 
comoving volume: 

ns = 5] vlnA . (2.2) 
A 

It is also useful to introduce the dimensionless parameter 

V = — , 2.3 
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which measures the baryon to photon number ratio. 

The mass fraction Xa of a given specie is instead defined as Q 



Xa = ^ (2.4) 



and verifies the normahzation condition 



Y.Xa = I . (2.5) 

A 



Finally, the abundance of a nuclide relative to hydrogen is 



Ya = ^ , (2.6) 
and, in terms of this, the mass fraction can be also expressed as 

In chemical equilibrium, the parameter fiA in (^3) is given by 

HA = Z lip + {A - Z) , (2.8) 

where fJ,p,fJ-n are the proton and neutron chemical potential, respectively. In terms of the 
quantities of the constituent nucleons, the number density ua can therefore be written as 



A-2 / 27r \ 5(^-1) 



B 



A 



p "n 



where mi\j is the nucleon mass and 

Ba = Z nip + [A — Z) rrin — ruA (2-10) 



is the binding energy of the given nuclear specie {Z, A). From ( |2.9| ) and ( 1.56| ) (applied to 



photons) we can then deduce the following formula for the mass fraction in terms of the 
nucleon quantities and the r] parameter: 

From (U) or (|2lll ) we then see that the greater is the binding energy of a given specie. 



the larger is its abundance (though, of course, this is not the only condition for having an 



appreciable amount of a nuclide). In table 2.1 we report the values of the binding energy 



of light nuclides relevant for nucleosynthesis. 

^In many papers, the mass fraction for the primordial ^ffe is usuaUy denoted with Y. 
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Table 2.1: Binding energies (in KeV) of light nuclides (data taken from |]35|| ). 



n 0.0 


76204.800 ± 0.421 


p 0.0 


^^C 73439.899 ± 0.952 


D 2224.573 ± 0.002 


79575.205 ± 1.400 


•i ZT 0/101 001 1 r\ r\r\ A 

H o4oi.ozi ± 0.0U4 


L 92161.753 ± 0.014 


'•^He 7718.058 ± 0.002 


'^'^N 74041.317 ± 1.000 


'^He 28295.673 ± 0.005 


^^C 97108.065 ± 0.016 


'^Lz 31994.564 ± 0.475 


^^A^ 94105.267 ± 0.270 


■^Li 39244.526 ± 0.473 


^''C 105284.507 ± 0.019 


■^Se 37600.358 ± 0.472 


^^A^ 104658.628 ± 0.016 


^Li 41277.328 ± 0.488 


^^O 98733.236 ± 0.076 


37737.813 ± 1.107 


^^iV 115491.930 ± 0.019 


30258.837 ± 62.471 


115955.627 ± 0.503 


64750.700 ± 0.370 


i^O 127619.336 ± 0.019 



2.2 Observed primordial abundances 

Primordial abundances of light elements are, in general, significantly altered by nuclear 
processing in stars; hence, the deduction of the primordial abundances from the ones 
observed today is very difficult and, in some sense, is always dependent by the chemical 
evolution model. 

Primordial D is easily destroyed by (p, 7) reactions in stars where the temperature is 
greater than 6x10^ it can be, furthermore, converted in ^He which is, however, 
also burned. Instead '^He grows in abundance with time: throughout stellar evolution, a 
star becomes more and more enriched in "^He and, through stellar mixing and mass loss, 
the interstellar gas can increase its content in ^He. Finally, '^Li is destroyed by nuclear 
reactions with protons in low mass stars, but it is nevertheless produced in supernova, 
red giant interiors, supermassive objects, etc. 

We will shortly review the actual experimental situation for the light element abundances. 
2.2.1 D 

Deuterium present in stellar interior is destroyed by the energy-generating reactions; in 
fact, it has not been detected in any star . Nevertheless, there are accurate data from 
giant planets in the solar system [^ ] and from the local interstellar medium (ISM) [^]. 
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Since there are no known astrophysical sources of D a firm lower bound to the 

primordial abundance from the lowest observed one can be obtained 

^ > 1.1x10"^ . (2.12) 

Other estimates involve models of galactic chemical evolution to infer the pre-galactic 
abundance in D from those actually present in the ISM. 

Deuterium has also been observed by many authors in the spectra of high red-shift quasar 
absorption systems (QAS) |]39[| , pO| , [41|. In principle these measurements would be 
capable of determining the primordial value oi D/H. However, at present, this is not the 
case. In fact, while in several measurements a rather high value of D/H ~ 2.0x10"^ has 
been reported ^M, for other observed QAS a significantly lower value (~ 3.4xl0~^) has 



been quoted |^, ||4l[. The present situation, then, does not lead to a conclusive value for 



D/H a "reasonable" bound from these observations is the following |42]: 

^ > 2.5x10-5 . (2.13) 
11 

2.2.2 ^He 

^He is the product of incomplete H burning in stars comparable in mass to the Sun, but 
it is burned away in the interiors of heavier stars. 

^He has been detected in meteoritic extractions [^] as well as in several interstellar 
medium (ISM) [Q| measurements of ^He in galactic HII regions Q and in planetary nebulae 
Observations show a wide dispersion which may be indicative of pollution or a bias 
, and large evolutionary uncertainties make very difficult to extrapolate the primordial 
abundance. However, noting that D is burnt in stars to ■^He a (fairly known) fraction of 
which survives stellar processing, a more reliable result can be obtained by considering 
the sum of D and "^He and requiring that these two isotopes not be overabundant at the 
time of formation of the solar system ||4^. This yields the limit 

^^^10- . (2.14) 
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Figure 2.1: The ^He (Y) vs oxygen (0/H) abundances in extragalactic HII regions (taken 
fromRef. 0). 



2.2.3 "^He 



The abundance of '^He has been determined in a variety of astrophysical sources |4S], [49|, 
such as atmospheres of young and old stars, planetary nebulae, HII regions (galactic and 
extragalactic) and so on. The almost, quite amazing constancy (within a ~ 20% factor) 
of the observed "^He abundance by mass of about 25 % in all these objects, as can be 



seen from Figure 2.1, points to an uniform origin, which is the primordial nucleosynthesis. 



However, although this abundance is measured with a much higher accuracy than for the 
other light elements, one has to be careful to extrapolate the primordial abundance from 
the today observed one because of stellar processing. In fact, stars generate their energy 
by burning hydrogen into helium and, at their death, they return this processed material 
to the ISM to be incorporated into subsequent generation of stars. 

To determine the primordial abundance of '^He, one must allow for the stellar helium 
component through its correlation with some other element which is made only in stars 
(such as N,0). The primordial value is then obtained with an extrapolation to zero 
metallicity. For this reason, it is better to concentrate on data from those regions least 
contaminated by the products of stellar evolution: the ^He abundance is, in fact, best 



■^These are regions in which interstellar gas is heated and ionized by the radiation from hot young stars. 
In astrophysical notations AI means neutral A element, All stands for A ionized once, and so on. 
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determined by observations of Hell Hel recombination lines in extragalactic HII regions. 
Measuring both "^He and 0/H abundances, it has been found that the data are well fitted 
by a linear correlation |5^] for Xi^^ (= Y) versus 0/H, and then the primordial value 
can be determined from the intercept of that relation. Moreover, the primordial "^He 
abundance can also be determined by studying the correlation between X^jj^, and N/H 
(in almost all HII regions both 0/H and N/H data are available). However, in this case, 
it is not still clear if a linear fit should be appropriate (at least in the situations in which 
N/H is not proportional to O/H) jS^, |52]. Nevertheless, the quoted difference 



between the intercepts of Xi^j^ versus N/H and Xa^^ versus O/H is small ( < 0.003). 
At present, there are two different quoted abundances |48|, [^] 

Xijj^ = 0.234±0.0054 (2.15) 
XiH^ = 0.243±0.003 . (2.16) 



The Izotov et al. |49] result (|2.5| ) is higher because of dropping the lowest metallicity galaxy 
from their sample due to its anomalously low Hel line intensities (which seems underlying 
stellar and interstellar absorption). Also, different ways to extract the primordial value of 
X4jjf, are employed. Although we do not enter into the details of the discussion for the 
quoted values, we limit ourselves to observe that the systematic error may be significantly 



larger than the ones included above [51|. In particular, the atomic physics inputs seem 
to be uncertain (regarding mainly the Hel line intensities) and more accurately known 
physical conditions in the HII regions are auspicable. This can be done by measuring 
several different line intensities. 

We finally note that the error on the "^He mass fraction is on the third significant digit, 
and consider the upper bound 

X,He < 0.24 , (2.17) 

as a reasonable one and 

X4He < 0.25 , (2.18) 

as a reliable bound [42|. 
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Figure 2.2: The ^Li abundance in halo stars with Tg^j-f > 5500 X as a function of metalUc- 
ity. Here [Li] = \og{^Li/H) + 12 and [Fe/H] = {\og{Fe/ H))/ {\og{Fe/ H)q) (taken form 
Ref. 0). 



2.2.4 ^Li 

As a star evolves, its surface lithium is subject to destruction and dilution. In fact Li, 
like D, is easily destroyed (at temperatures above ~ 2x10^ K) by (p, a) reactions, so that 
only the lithium remaining on the stellar surface survives. This, however, will be further 
diluted through mixing of the outer layers with the interior. To avoid this, stars with 
a surface temperature T > 5500 ii' and a metallicity less than l/20th solar are mainly 
observed, so that effects such as stellar convection may be not important. 
"^Li is observed in the atmospheres of both very old stars (Population II) and young ones 
(Population I) [ |53| , p^]. When we plot the ^Li abundance as a function of metallicity 
(for stars with T^^jr/ > 5500 K) a plateau region is observed, as can be seen in Figure \2.2[ 
indicating a primordial (or very close to it) value for the abundance. 
The best estimate for the mean ^Li abundance (including statistical uncertainties) is the 
following 

(1.6 ±0.1) X 10"^° . (2.19) 



H 



However, we stress that Li abundance determination is sensitive to several stellar pa- 
rameters (surface temperature, metallicity and so on) and an important source of sys- 
tematic error is due to the possibility that Li has been depleted in stars from their 
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initial abundances. These uncertainties are limited by the observation of ^Li [55| (with 
^LiP Li = 0.05ib0.02), a very fragile isotope; in fact standard stellar models predict that 
any depletion of Li would be accompanied by a very severe depletion of ^Li [p^. 
Finally, Li is also produced, along with Be and B, in cosmic ray spallation of C, N ^ O 
by p, a (and also by a — a fusion). Hence abundances of Be and B have been served as 



a consistency check |57] on primordial Li, concluding that no more than 10 - 20 % of the 
Li is due to cosmic ray nucleosynthesis. 

Of course, for Li the uncertainties are dominated by systematic effects. 



Abundances of other elements have also been observed, such as those of the inter- 
mediate mass isotopes ^Se, ^^B and ^^B. However, for these, large uncertainties are 
introduced, since it is believed that these isotopes are formed in cosmic ray nucleosynthe- 
sis. In fact, the observed abundances are far above the BBN predictions and a comparison 
between theory and experiments is extremely difficult. 

2.3 Production of the light elements 

Before going into the details of the primordial nucleosynthesis, in this section we will 
discuss (mainly qualitatively) the main features and results of Big Bang Nucleosynthesis. 

2.3.1 T>lMeV{t<ls) 

At sufficiently high temperatures, neutrons and protons are maintained both in kinetic ^ 
and chemical equilibrium by charged current weak interactions: 

Ve + n ^ e~ + p (2.20) 
e+ + n ^ p + T7e (2.21) 
n ^ p + e~ + Ve . (2.22) 

Because of chemical equilibrium, the following relation between the chemical potentials 
holds 

Hn + t^u = IJ-p + fJ-e ■ (2.23) 

■^Kinetic equilibrium is maintained also and primarily by nuclear and electromagnetic elastic scattering 
reactions 
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The neutron to proton density ratio is then given by 



= e T e T (2.24) 



Up I 

. ^ / eq 



(detailed balance), where Q = rUn — rup ~ 1.29 MeV is the difference between the neutron 



and proton mass. As discussed in sect. |l.3.l| , the electron chemical potential is small 



compared to the temperature T in the region relevant for BBN. Furthermore, assuming 
zero neutrino chemical potential, we have then 



The nucleons abundances track their values in equilibrium, 



Xn{T) ^ Xp{T) = , (2.26) 

eT +1 



as long as the rates for the reactions ( 2.20| ),( p^ ),(2.22) decrease sufficiently and become 



comparable with the expansion rate H. This happens, as we have seen in sect. 1.4, at 
about T ~ 1 MeV. At this time, neutrons "freeze out", i.e. go out of chemical equilibrium, 
and Xfi (and Xp) relaxes to the final constant value 

X^ ^ , (2.27) 

e'^F + 1 

(Tp ~ 1 MeV is the freeze out temperature) rather than following the exponential falling 
of ( 2.261 ). Since Q/Tp ~ 0(1), a substantial fraction of nucleons survives when chemical 



equilibrium between them is broken, while if Tp would be much lower than Q then n„ ~ 0. 
The abundances of the other nuclides are very small for T ^ 1 MeV, and can be calculated 
from ( 2.111) . They become appreciable when the temperature goes down the binding energy 



per nucleon of the given specie, which is typically of the order of 1 8 MeV. However, 
there are two facts that delay the onset of nucleosynthesis: the low binding energy of the 
first nuclide in the nucleosynthesis chain, deuterium, and the high number of photons (or 
high entropy) present at the nucleosynthesis era. In fact, for T < 1 8 MeV deuterium 
(and hence the other nuclides to be formed from this) synthesis is energetically favoured, 
so the reaction n + p D + j takes place but, due to the high number of photons, 
the inverse reaction takes place as well and more efficiently than other reactions such as 
n + D ^ + 7 or p + D ^ ^He + 7. Hence the abundances of these light elements are 
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still very small for these temperatures (deuterium "bottleneck"). For reference, we report 
on the typical abundances of some light nuclei for T ~ lOMeV, as calculated from ( 2.11| ) 
with r] ~ 10~^: 

X„ ~ ~ 0.5 , (2.28) 

Xd ^ 6x10-12 ^ (2.29) 

Xbh^ ~ 2x10-23 ^ (2.30) 

XiHe - 2x10-3^ , (2.31) 

Xi2c ^ 2x10-126 . (2.32) 



Instead, as neutrons freeze out, the n/p ratio is given by (2.25) with T ~ Tp {rin/np ~ 1/6) 
and for T ~ 1 MeV the abundances of light nuclides are 



Xn 




1 

7 ' 


(2.33) 


Xp 




6 

7 ' 


(2.34) 


Xd 




10-12 , 


(2.35) 






10-23 , 


(2.36) 






10-2*^ , 


(2.37) 


Xl2c 




10-108 _ 


(2.38) 



2.3.2 r~0.3^0.1Mey (t~ 1^3mm) 

As the temperature cools below 0.5 0.3 MeF the synthesis of complex nuclei becomes 
thermodynamically favourable. Thus, neutrons and protons react with each other to build 
up light nuclides through the following sequence of two-body reactions: 
p{n,-/)D 

D{p,'y)^He D{D,n)^He D{D,p)^H 
^H{D,n)'^He {'^He,-/)'^ Li 

^He{n,p)^H ^He{D,p)^He ^He{'^He,-f)'^Be 
^Li (p^He) ^He ''Be {n,p) Li 



(2.39) 



(a complete list of nuclear reactions is reported in appendix |^). The first reaction is, of 
course, crucial since deuterium has to be formed in appreciable amount, before the other 
reactions can start. In fact, many body interactions such as 2p + 2n -^'^ He are in general 
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ineffective, due to the small number densities of nucleons. After deuterium, appreciable 
quantities of '^H and "^He may form. The synthesis of ^i7e is, however, delayed until 
enough tritium has been built up, since the main process for making '^He involves this 
hydrogen isotope. 

Nearly all neutrons surviving down to freeze out are captured in '^He because of its large 

binding energy per nucleon; its mass fraction is then approximately given by 

4n4He _ 4(nn/2) 2{nn/np) 
-^^He = - \ = -rn 1 — r ■ (2.4U) 

For T ~ 0.1 MeV, the n/p ratio is about 1/7 (actually, it is a bit smaller than the freeze 
out value since neutrons have been depleted by /3-decay), so that 

X,He ^ \ ■ (2.41) 

Heavier nuclei do not form in any significant quantity both because of absence of sta- 
ble nuclei with ^ = 5,6 which forbids nucleosynthesis through {n,^ He), {p,^ He) or 
{^He,'^ He) reactions) and the large Coulomb barriers for reactions such as ^H(^He, Li 
and ^He{^He,-f)^Be. 

For t > 10^ s BBN is over, since the low temperature and low density in the Universe 
suppress nuclear reactions. Hence the most abundant nuclides with which we are left are 
hydrogen (p) and '^He (incorporating almost all neutrons), followed by trace amounts of 
D, ^He and ^Li. Note that the "^He abundance includes that of survived ^H which sub- 
sequently /3-decays and similarly the "^Li abundance includes that of ^ Be. 
Observe that the ^He abundance, as approximately given by ( p. 41 ), is quite insensitive to 
the nucleon density (or rj) present in the Universe, and almost only depends on the weak 
reaction rates determining neutron freeze out. This is not true for the other elements, 
for which the nucleon density directly determines the two-body nuclear reaction rates. 
It is expected that D and ^He abundances decrease with increasing ry, since for higher 
values of tj we have a more efficient burning into '^He. In fact, from ( p. 11 ), abundances 
increase with r] so that BBN can start earlier when at higher temperature and the n/p 
ratio is larger. As far as "^Li abundance, it is expected to decrease with increasing rj when 
it is determined by the competition between '^Hei^ H,^y Li and Li{p,'^ He)^He, while 
for sufficiently high rj it starts increasing again with r] due to larger production of ^Be 
through ^He{^He,^y Be{e~ ,Vey Li. 
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Figure 2.4: The intermediate mass element abundances from big bang nucleosynthesis as 
a function of r] (taken from Ref . |^8| ) . 



In Figures 2.4 we report ||T^ the predicted primordial abundances of some light ele- 
ments versus the baryon to photon ratio rj as calculated with the standard numerical code 
[p9|]. The discussed features come evident from these plots. 



2.4 Calculation of the primordial abundances 

In this section we will write down the equations governing the time evolution, during 
nucleosynthesis, of the abundances of the produced light elements. 

The variables describing the Universe as a thermodynamic system during BBN are the 
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following: 

R 

T T 

(2.42) 

Yi . 

The first is associated to the dynamical aspects of the system, i.e. expansion of the Uni- 
verse, while the others describe the primordial plasma of e±, i/, n,p and Yi are the nuclide 
abundances. 

The relevant part of nucleosynthesis takes place after neutrino decoupling, so that the 
quantities describing the neutrino plasma evolve independently from the plasma in equi- 
librium and, for our purposes, may be disregarded. Hence, Ty and fiy may be dropped 
from the list above. 

Moreover, only two of the remaining three chemical potentials are independent, since the 
relation ( p. 23 ) holds at chemical equilibrium. The first independent variable we choose is 



the electron chemical potential or, equivalently, the parameter (p^, = /Ug/T. Instead, the 
second one, which is related to the nucleon chemical potential, is the total baryon number 
density ns- Thus the independent variables we have to consider are the following 

i?, ns, (/.e, T, , (2.43) 

for which it is necessary to write down the time evolution equations; however, let us note 
that in many cases is more convenient to follow the temperature evolution rather the time 



one, and this can be done once the relation between T and t is established |60, 61, p9| . 
2.4.1 Equation for R 

As seen in the previous chapter, the equation satisfied by the cosmic scale factor is the 
following: 

where p is the total energy density. During the nucleosynthesis era, the curvature term 

k/R'^ is completely negligible, so that we can safely neglect its contribution to expansion. 

*Note that it is explicitly assumed that species interacting electromagnetically or strongly are always 
in kinetic equilibrium, though chemical equilibrium may be not hold. For neutrinos, it is instead supposed 
that their decoupling is "instantaneous" (see below and chapter H) 
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2.4.2 Equation for ub 

The total number Nb of baryons in the comoving volume is constant ^ and, assuming 
homogeneity, it is also independent of R. Hence the baryon number density ub = Nb/V 
varies as , so that 

Ub at K at 

which is the equation for the evolution of n^. 

2.4.3 Equation for 0e 

The electron chemical potential (or (j)e) is fixed by the conservation of the electric charge 
and the total neutrality of the Universe [pOl 



ne- - = Up , (2.46) 
neglecting the possible presence of antiprotons. The proton number density is 



Hp = ^ ZiTli = UB Zi Yi = UBQB , (2.47) 

i i 

while, for 4>e <^ 1, we have 



ne- - n,+ ~ l^ct>eT'^ f{z) , (2.48) 



f{z) ^ / dx 



oo ^2 







(e^ + 1)2 



with e = y/ x"^ + z"^ and z = me/T. From these we then obtain the equation for 



~ (2.49) 



2 T^f{z 



with qB = Y.i ZiYi defined in (|2^ . 
2.4.4 Equation for T 

The important equation relating temperature with time is deduced from energy conserva- 
tion, which gives 

|{.V)+p^ = , (2.50) 



^At the energy scales relevant for BBN there are no baryon number violating processes. 
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where V ^ is the comoving volume and p, p are the total energy density and pressure 
of the particles involved, namely: 



P = P-Y + Pe + PB 
p = P7 + Pe + PB 



(2.51) 
(2.52) 



(assuming that neutrinos are completely decoupled ^ The photon contributions P"f,P"f 



depend only on the temperature and are given by equations ( 1.57 ), ( 1.58| ), while pe 
Pe- + Pe+1 Pe = Pe" + Pe+ depend in general on both T and (pe'- 



Pi 



vr 
15 



1 



P7 



e 



giz) = I dx 
^ ' Jo + 1 



h{z) 



Pi 



dx 

e e^ + 1 



(2.53) 
(2.54) 
(2.55) 



The baryon terms can be written, in the non relativistic approximation, as (see ( |1.52| ), 

(S)) 

PB = Pi ('^i'^i ^ ^"^""V ^ 

i i 



PS 



Mu 2 Mu. 



(2.56) 
(2.57) 



where M„ is the atomic mass unit (referred to ^^C) while Amj is the mass excess for the 
nuclide i (mass excesses for the light nuclides relevant for BBN are reported in appendix 



From (2.50) after some simple algebra we arrive at the formal equation for T 



dT 
Itt 



I dR (dp-'^ 



(2.58) 



First of all, let us note that the energy density derivative with respect to T has to be 
evaluated along the time direction, that is 



dp _ /dT 



\ 1 dp (dT\ ^ fdp^ dpe dpB\ , . 



^In this approximation, neutrino plasma energy is separately conserved, 
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dp^ dp^ dT 

dpe _ fdpe dp^d(j)e\ cff / dpe d(j)e duB dpe d(j)e dqs 

dt ~ \dT d(l)edT) dt \d(l)ednB dt d(l)edqB dt 

dpB (_o]_^ , 3 J_dr ^ , V- V. ^ m ^ 

dt ^ ^^1^ '^Rdt + 2 M„ Y ^ i V 2MuJ dt 



Hence we have 

dp 1 dpe dpe d<j)e dpj 3 1 

dT ~ \~df ^ W^~dT ^ It ^ 2M, 



ldRfdT\-^f ( ( IdRy^^fAmi 3 T \ dVA 

^ 1 dR\^^ / dpe dcpe dns ^ dpe dcpedqB 



R dt J \d4>ednB dt dcp^dgB dt 
and substituting Eq. ( |2.58D , solving with respect to dT/dt, we find 
dT r„ 1 dR 



dt 



( 1 dR , 

I 3 — — (p^, + Pe + P7 + Pe + Pb) + 



dpe dcpe duB dp^d(i)^dqB_ ^ f Am^ dYj \ 

d(t)ednB dt d(t)edqB dt ^ \ Mu 2M„; dt J 

J dpe dped(l)e dp^ 3 I ST v\ ^ fo «n^ 

This is the equation relating temperature with time we looked for. Note that it is now 
written in terms of known quantities, once Eqs. dHI) and are taken into account. 

2.4.5 Equations for Yi 

The change in time of the abundance Yi of a given nuclide i is driven by the rates for the 
nucleon reactions producing or destroying it which we will consider to be of the form 

i + j < — > k + I . (2.61) 

If Tij_ki is the reaction rate for the direct process and T^i^ij that for the inverse one, the 
time evolution of Yi is given by the following equation: 

^ = Tki^,,YiYk - Tij^klY^Yj . (2.62) 

In general, one can consider also more complex reactions, in which more than one nuclide 
for each specie is present: 

Ni (i) + Nj (j) ^ Nk (k) + Ni (0 , (2.63) 
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Figure 2.5: '^He mass fraction and D, Li abundances predicted in the standard BBN 
scenario. Also shown are the observed values of these isotopes (adapted from Ref. 0]). 

where Na is the (integer) number of a given nuclide taking part to the considered reaction. 
In this case the evolution equations are of the following general form: 



dYi ( y^iy^fc y.^» y.^J 



(2.64) 



The complete nuclear reaction network used in BBN calculations is reported in appendix 

0- 

2.5 Theory versus observations. Constraints from BBN 

We can now compare the inferred bounds on the abundances of light elements discussed 



in section "L2 with their values computed within the standard BBN scenario. 

The two conflicting '^He mass fraction results reflect into the QSO determinations of 

deuterium abundance, and two different sets of measurements mutually incompatible 

XiHe = 0.234±0.0054 , D/H = (1.9±0.4)-10-^ (2.65) 
XAHe = 0.243±0.003 , D/H = (3.40±0.25)-10"^ (2.66) 
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arise. Obviously, large helium mass fractions require low values for deuterium abundance 



and vice versa. The situation is well illustrated in Figure 2.5 where both theoretical 
predictions and experimental observations are reported; from this the ranges of the r] 
parameter corresponding to the different sets of data (2.65), (p.66|) can be directly read 



off. However, a very recent (re-)analysis of observational data has superseded the 
discussed dichotomy, and a unique range for r] emerges, as can be seen from Figure |l[ 
Nevertheless, adopting the reliable bounds of Eqs. ( 2.13| ), ( ^.181) , (|2.19| ) we can derive the 
conservative limit on r/ 

4.1x10"^^ < 1] < 9.1x10"^° , (2.67) 

overwhelming all the quoted difficulties in abundance determinations. These values of 
r] translate into the following range for the nucleon density (normalized to the critical 
density) 

0.0015 < nsh"^ < 0.033 , (2.68) 
which is then an independent estimate of the baryon content of the Universe. Confronting 



with Eq. ( |1.36 ), it seems evident that non baryonic (dark) matter exists in the Universe. 
For an account of dark matter problems, see for example []63| . 

Having established the consistency of standard BBN, we can use it to constrain new 
physics beyond the Standard Model of elementary particles. Here we don't perform an 
exhaustive analysis of these constraints, which is far from the scope of this thesis and for 
which we remind to excellent reviews (such as Q); our aim is to show how precise BBN 
results can be used to extract important informations on particle physics. 
Limits on this are mostly sensistive to the bounds imposed on the ^ife abundance, which is 
predominantly determined by the neutron to proton ratio at freeze out. This is determined 
by the competition between the weak interaction rates and the expansion rate of the 
Universe, according to Eq. (|1.107| ). The presence of additional neutrino flavours (or any 
other relativistic particle species) at the time of BBN increases the overall energy density 
and hence the expansion rate, leading to a larger value of Tp, n/p and ultimately X4fff,. 
The dependence of Xajj^, on the effective number Ni, of neutrinos can be parametrized as 
follows @: 

X4He = 0.2262 + 0.0131 {N^ - 3) + 0.0135 In (j^^ (2-69) 
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(this fit holds near the central value of ( 2.67 )). Note the weak logarithmic dependence 
of XiUf, on the r] parameter. By requiring that all constraints on elemental yields be 
simultaneously satisfied, a Monte Carlo analysis has given the following conservative limit 
on ig 

< 3.75 + 78 {X^He - 0.24) . (2.70) 

A number of other constraints on known or hypothetical particles, both relativistic 
and non relativistic at the time of BBN, as well as limits on the strength of new forces, 
can be deduced. These also have important implications for the problem of dark matter 
in the Universe. While these constraints are much more stringent than the ones obtained 
in laboratory experiments, pointing out the relevance of BBN studies, they are (particle 
physics) model dependent and cannot be discussed here. We remind to the recent reviews 
[ p^ , |7| and references therein. 
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Chapter 3 

The Born rates for n ^ p reactions 



The initial conditions for the primordial nucleosynthesis are settled out by the ratio of 
the relative abundances of neutrons and protons. When BBN is starting, the primordial 
plasma consists of nucleons, electrons, positrons, neutrinos, antineutrinos and photons in 
thermal equilibrium. The relative abundance of neutrons and protons is determined by 
the following charged-current weak interactions: 

(a) Ue + n ^ e~ + p , (6) e~ + p ^ Ue + n (3.1) 

(c) e'^ + n ^ Ve + p , (d) Ve + p ^ e'^ + n (3.2) 

(e) n ^ p + e~ + Ve , {f)p + e~+Ve^n (3.3) 

In this chapter, the leading contribution for the rates (per nucleon and per time) for these 
six reactions is calculated and the relevant approximations are particularly pointed out. 
For a generic process in ( |3.1| ), ( p. 21 ), ( p.3| ) , the reaction rate can be written in the form 

(3.4) 



spins \ ct / 



spins 

|2 



where |M| is the squared matrix element, which has to be summed over all spin degrees of 
freedom. The 5-function gives the energy-momentum conservation, p^ = {Ea, Pa.) being 
the 4-momenta of the particles entering in the process (each pa is intended with a positive 
sign if the particles involved in the reaction is in and with a negative sign if it is out). 
The quantity = Ha ^^^a is the phase space factor for the given reaction. Since the 
processes (|3.lD, (p. 2D, (3.3) occur in a thermal bath, the phase space factor associated with 



each initial (Fermi) particle is 



= ^^J_Fa , (3.5) 

(27r)3 2Ea 
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while for final (Fermi) particles we have 

where Fa is the phase-space density of the a-th particle at temperature T, which takes the 
thermal equilibrium value (in the rest frame of the thermal radiation) given by the Fermi 
function 

FaiEa) = ^ (3.7) 



only for those species actually in equilibrium. Note that in writing p.4| ) we have assumed 
that Ua does not depend on spin, so that the summation over this degree of freedom is 
limited to |Mp. 

The quantity of interest in BBN is the reaction rate per incident nucleon, so that F in 
( |3.4D has to be divided by the number density 

= ^ y (# ^^-^^ 

of the incident nucleon (say b). Assuming that X^spins \M\'^ does not depend on the ini- 
tial nucleon momentum (this assumption, which is justified in the infinite nucleon mass 
approximation, will be discussed later), the net effect on the reaction rate is the following 

dm - /) = ^ E (2^)' (e P-) n ^ (3-9) 

spins \ a / aj^b " 

(from now on, we indicate with F the reaction rate per incident nucleon). Furthermore, 
since at the epoch of BBN the baryon density is very low (r/ ~ 10~^), we can also neglect 
the occupation number of the final nucleon (say c), so the final expression for F(i f) to 
be considered is 



dr(i^f) = l Y. \Mf i2.)U^ W n ^ ^ ■ (3.10) 

spins \ a / \j,eptons / 

We now proceed to explicitly calculate the reaction rates for the processes in ( |3.1| ), ( |3.2| ), 
(|3.3D in the Born approximation. 
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Figure 3.1: The Feynman diagram for the reaction + n ^ e + p. 

3.1 The reaction + n ^ e~ + p and the other processes 



Let us start with considering the direct reaction in ( |3.lD ; in the Born approximation the 
Feynman diagram for this process is given in Figure The amplitude is 



Gp _ 



= — Up(P4)7m(Cv - C'A75)^in(P2)^ie(P3)7''(l - 75)^ii/(Pl) , (3.11) 



where Gj? is the Fermi couphng constant and Cy, Ca are the vector and axial coupling of 
the nucleon. In the neutron rest frame (which, in the limit ttt-at — > oo, coincides with the 
plasma rest frame), neglecting proton motion (p4 ~ 0, £'4 ~ -/W4), we then have 

|M|2 = 32G| M2M4 [(Cy'' + ?,Ca^) EiE-i + (c^' - Ca') P2Pz cos 6,,) (3.12) 

spins 

(hereafter pa = |Pa|)i In r(i/e + n — > e~ +p), after integration over angles, the last term 
does not contribute, thus we can disregard it, and obtain the well known result 

^ |M|2 = 2,2Gl (C7^2 + 3C7/) M2MiEiE^ . (3.13) 

spins 



The reaction rate in ( 3.10| ) for the present process then takes the form 
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■ {2Tr)^ 5^ {pi + P2 - P3 - P4 



3 • 

1111 



2M2 2M4 2Ei 2E3 

d^pi d3p3 c^Sp^ , . , . 

(we have distinguished the neutrino Fermi function from that of electrons since, in general, 
neutrino temperature Tj, is different from the 6^,7 bath temperature T, as discussed in 
the previous chapter). The integration over the proton 3-momentum can be eliminated 
by using 5'^(pi + P2 — Ps — P4)- From the isotropy of space, also the integration over the 
incident neutrino solid angle / di^i^^ can be easily performed, leading to an overall factor 
47r: 

dT{ue +n ^ e-+p) = ^ ^ 5{Ei-Q)pjdpipldp3 dfl^ F^{pi) (1 - Feips)) , 

(3.15) 

where Q = M4 — M2 + E3. The energy (5- function can be used to perform the integration 
in dpi, assuming massless (or almost massless) neutrinos, i.e. Ei ~ |pi|- Finally, also the 
integration over the electron angles / df^s is easily performed since the differential rate 
does not depend (in the present approximations) on these quantities. We are then left 
with the following result; 



r{ue+n^e-+p) = ^^^^3 ^ J dp3plQ^9{Q)F,{Q){l-FM) . (3.16) 

The integration limits are imposed by the condition Q > 0; in the present case we have 
P3 G \/ — m^, +co^ , with A = M2 — M4 (in the following, with A we always denote 
the neutron - proton mass difference). 

For the other five processes in ( |3.lD , ( |3.2| ) , ( |3.3[) one simply has to observe that the form 



of X^spins \M\ in (3.13) remains unchanged due to crossing symmetry. Then in the final 
result (3.16) we have just to replace Q (which is determined by the energy conservation for 
each reaction) and the thermal factors (which depend on the initial and final lepton states 
for each process). In Table |3.l| we report the values for Q and the thermal factors for each 
reaction in (|3.l|) , (^), (|3.3D , using the same notation adopted for i^e + n — > e~ +p (i.e. 
1 = Ue or I7e , 2 = n, 3 = e~ or e"*", 4 = p). The Born rates for each reaction in ( |3.lD , (|3.2| ), 
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Table 3.1: The Q parameter and thermal factors for the reactions in ( |3.lD , ( |3.2[) , (3.3) 



Reactions 


Q 


Thermal factors 


Uf. n ^ e~ p 
e~ p — > Vf, n 
n — > I7e p 
Ve, p — > n 
n pe~Ve 
pe~Ve — > n 


-A + E3 
-A + E3 
A + Es 
A + E3 
A - E3 
A - E3 


F,{Q) (1 - F,{E3)) 
FeiEs) il-F,iQ)) 
FeiEs) (l-F.iQ)) 

F,{Q) (1 - Fe(i?3)) 

(1 - F,{Q)) (1 - FeiE^)) 
F,{Q) FeiEs) 



(3.3) , as functions of photon temperature in the range relevant for nucleosynthesis, are 



plotted in Figure 3.2. The relation between neutrino and photon temperature we adopt is 



reported in the Appendix ^ 

3.2 Radiative electromagnetic corrections: the neutron life- 
time 



The reaction rates for the six weak processes in ( p.ip , ( |3.2[ ), (3.3) are proportional to the 
quantity 



Gl (C7^2 ^ 3^^2^ =Gl. [1 + 3^ ) , (3.17) 



as it can be seen from Eq. ( 3.16| ), for example. The value of Gp is very well known from 



the measurements of the muon decay rate 111] , while Gy is obtained from the observations 



of 0'*' 0'^ nuclear (3 decays after applying radiative and isospin-mixing corrections [18|. 
Instead, the empirical value of the ratio Ca/Cv is deduced, for example, from experiments 
measuring the angular distribution of decay products of polarized neutrons [nsl. It is very 



common in the literature on BBN (see for example Q) to rewrite the effective coupling 



constant in (3.17) in terms of the neutron lifetime (free neutron decay in vacuum) since. 



as it can be immediately deduced from ( 3.16| ) and Table ^T] dropping out the thermal 



factors, it is given by the following relation 



^'"^ dee (^e - (^^ - l) ^ , (3.18) 



where mg is the electron mass. 

Rewriting the reaction rates in terms of the neutron lifetime, it is then customary to 
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T/m. 



Figure 3.2: The Born rates Tb versus photon temperature divided by electron mass nie- 
Hereafter reaction labels are the same as in ( |3.2D , ( |3.3D . 



substitute with its experimental value. This way of reasoning does not bring to accu- 
rate BBN theoretical predictions at a certain degree of approximation, but introduces an 
uncertainty in the BBN calculations due to the experimental uncertainty on the neutron 
lifetime. This procedure is usually adopted to circumvent the problem of an accurate 
estimation of radiative (electromagnetic) corrections to Gy which is, nevertheless, a very 
important one. To give an idea of the relevance of this problem, it suffices to observe 



that if we calculate the neutron lifetime from the tree level formula ( 3.18| ) (not containing 
radiative corrections), inserting the value of Gy and Ca/Cv reported in [|l8|, we obtain for 
Tn the value ~ 961 s which is not at all consistent with the experimental value 886.7 it 1.9 s: 
we need a correction factor of about 8 %. 

This simple exercise shows that the inclusion of radiative corrections is a very important 
step, which we are going to discuss. 

Let us consider order a corrections to the neutron lifetime (a is the fine structure con- 
stant). These can be separated into "outer" corrections, involving the nucleon as a whole, 
and "inner" corrections, depending on nucleon structure. Obviously, inner corrections are 
sensible to the details of the strong interactions inside the nucleon, while it can be shown 
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[|6^ that outer corrections, at least up to terms of order a, are independent of these. 
The main Feynman diagrams contributing to order a corrections to neutron Ufetime are 



sketched in Figure 3.3. The electromagnetic interactions of the proton are through its 
charge and magnetic moment, while those of the neutron are only through its magnetic 
moment. All diagrams reported in Figure are of order a but the ones involving mag- 
netic moment interactions of the nucleons are suppressed with respect to the others by the 
inverse of nucleon mass. Notice that since we use the value of Gi? as measured from muon 
decay, we have also to consider radiative corrections for the muon decay rate (at the same 
order in a). The relevant diagrams are similar to those pictured in Figure |3.3| , but with 
n,p replaced by /i, ; obviously for muon decay there are only outer corrections, since we 
deal with pointlike particles. In the following we consistently take into account also order 
a corrections to muon decay and report the global correction factor for the neutron decay 
rate. 

The outer correction to the (non integrated) decay rate can be written as 

(X 

— g{E,E^) , (3.19) 

where g{E, Em) is a function of electron energy E and end-point energy E^m reported in 
This function, which is the same for electron and positron capture, describes the 
deviations from the allowed electron spectrum arising from the radiative corrections of 
order a |^. It increases the decay probability for neutron /5 decay (and then decreases 
neutron lifetime) by about 1.5 %. 

Inner corrections are much more difficult to handle, since they strongly depend on nucleon 
structure. In general, one can follow two different approaches. On one side one can directly 
consider radiative corrections to the effective nucleon weak current Up7^(l — ^75)un | p5[ 
while on the other side, one can study corrections for the weak quark current ^7^(1 — 75)^ 
and then translate the quark-based description to hadronic description |^6| . Here we adopt 
the second point of view, and report the results obtained by Marciano and Sirlin [pGl 



^It is common to substitute the function g{E,Em) in (S.19) with its mean value 'g{Em) obtained 
averaging over the allowed elect ron spectru m. H owever, since we will consider not only neutron decay but 
also the other five reactions in (|3.l[), (p^), (3.3) , we do not employ this approximation. 
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The first term represent the dominant model-independent short-distance contribution (Mp 
is the proton mass, while Mz is the Z boson mass). The second and third terms are 
axial-current induced contributions, where is a low energy cutoff applied to the short- 
distance part of the box diagram and 2C is the remaining long-distance (low energy) 
correction. These terms depend on the details of the strong interaction structure and are 
the main sources of uncertainty in the radiative corrections. The allowed range for the 
cutoff Ma is 400 ^ 1600 MeF, while C = 3(7^ •0.266- (/ip + fin), fJ'p + IJ'n ^ 0.88 being the 
nucleon isoscalar magnetic moment. We then have 

— (ln^ + 2C) ~ 0.0012 ±0.0018 . (3.21) 



27r V Ma 



The last term in (3.20) is a perturbative QCD correction whose calculation is rather 
reliable and gives 

Ag ~ -0.34 . (3.22) 

Let us observe that the largest correction comes from the first term in ( 3.2C| ). Thus, it 
seems appropriate to approximate the effects of higher orders by summing all leading- 
logarithmic corrections of the type a" In" Mz {n = 1, 2, 3, ...) via a renormalization group 



analysis. This has been done in [37|; the corrected rate then acquires the factor 
Q{E) = (l + ^ (in^ + 2C) + {9{E,Em) + Ag)^ S{Mp,Mz) . (3.23) 



where a{fi) is the QED running coupling constant defined in the MS scheme satisfying 
the equation 

fi^aifi) ~ boifi)a\f,) , (3.24) 

with 

^o(m) = E ^3/ - ^f) - ^ - • (3-25) 

The sum runs over all elementary fermions with mass mj and charge Qj, Mw is the W 
boson mass. The short-distance enhancement factor S{Mp, Mz) is given by [^], |66] 



3 



S{Mp,Mz) 



JL , . JL , , JL , . 36 

17 



a{mc) \ * f a{mr) f a{mb) a{Mw) f a{Mz 



a{Mp) I \a{mc) / \a{mr) ) \ a{mh) J \a{Mw 



(3.26) 



^The expression reported in is here corrected allowing the (now observed) top quark mass to be 
higher than the W, Z boson masses 
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Table 3.2: QED running coupling constant for several energy scales. 
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In the MS scheme one has a~^{Mz) — 127.90 [|T8|. The running coupling constants at 
different scales are reported in Table |3.2| where we use the central values for the quark 
masses as quoted in [|l^]. From this one easily finds 

S{Mp,Mz) ~ 1.02254 . (3.27) 



The radiative corrections considered here, which are described by the factor in ( 3.23| ), 



are not, however, the only corrections to be taken into account for the neutron decay (they 
contribute for about 4 % to the lifetime, against the required global correction of about 
8 %). The remaining leading corrections are usually viewed as corrections to the phase 



space factor namely to the integrand function in ( 3.18 ) since they are (electron) energy 
dependent. We will express the contribution to the phase space factor as multiplicative 
(energy dependent) terms to the integrand function. 

Let us first examine the correction due to the distortion of the outgoing electron wave 
by the Coulomb field of the proton. This can be calculated by solving the Dirac equation 
for an electron under the influence of a spatially finite proton charge distribution of radius 
i? ~ 1 fni; the wave function is then evaluated at the centre of the proton. The correction 
factor is then given by |68| 

HE)CiE) + an^=^=^) (l " + l)) • (3-28) 



Note that the first term in ( |3.28| ) is usually denoted as the Fermi function for the Coulomb 
scattering. 

Another relevant correction comes from the fact that neither the electron wavefunction 
evaluated at the centre of the proton through T{E)C{E)) nor the (anti-)neutrino one 
are constant through the proton volume. Thus the decay rate has to be found by an 
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appropriate convolution of the electron, (anti-)neutrino and proton wavefunctions through 



the proton volume. This brings to the following correction factor [68 



C ~ 1 + Co + ^ + C2 ^ + C3 ^2 , (3.29) 
hi 



Co = 




3 


Cl = 


2i? „o ^ 
— ml Era 
15 




C2 = 


^(3 - B)R 


2 


C3 = 


A(3-i?)i?2 




B = 







Q 



aBREr, 



[Cy^ + 3C^ 

Finally, there are also small corrections related to the fact that the outgoing proton is 
not at rest, since it has a finite mass. This type of effect will be considered in detail in 
the following chapter. Here, we only want to point out that if the proton recoils, then the 
Coulomb field that distorts the electron wave comes from a moving source. The correction 



associated to this effect is contained in the factor [S8| 



From previous results, the neutron lifetime corrected at order a against electromagnetic 
interactions reads 



G'^(Cy^ + 3CA rA 2 / , 

^-fdEEiE-Af [E^-mt)^ 



g{E)J^{E)C{E)C{E)Q{E) , (3.31) 



where Q{E) is the correction function in ( p. 23 ) and the other terms are reported in ( 3.2g| ) 



- ( |3.3C1| ). Concerning the relevance of the different corrections just considered, a comment 
is in turn. For seek of completeness we have reported an exhaustive description of all 
corrections at order a, but as clearly appears from the explicit expressions of ( |3.23 )-( p.30 ). 



the main contributions come from ( |3.23| ) and J'{E) term in ( |3.28| ). All other contributions, 
contained in C{E), C{E) and Q{E), are in fact much smaller and thus they can be safely 
neglected, since they are suppressed by a factors of the order A/Aqcd- Actually the 
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contributions coming from radiative and Coulomb effects are even larger, or of the same 
order of magnitude, than these terms, so they should be included for consistency if C{E), 
C{E) and Q{E) are taken into account ||6^. For the level of accuracy of our analysis it 
will be sufficient to include J~{E) and Q{E) at order a only. 



Evaluating numerically the integral in (3.31) we then obtain for the neutron lifetime the 



value 893.8,5 which is now quite compatible with the experimental value. Note that the 
theoretical prediction can be further refined (leading to a slightly better agreement with 
the experiments) by considering other small effects (again of order a) such as magnetic 
moment (normal and anomalous) interactions, residual average proton polarization due 
to parity non-conservation and so on. While these effects are briefly discussed in [^^, for 
our purposes we do not consider them, since their contributions are energy independent 
and then account only for an overall factor. 

3.3 Results for the reaction rates 

The analysis of radiative corrections for /3-decay shows that their contribution to the 
rates for the processes in (0), (|3.2|), ( ^.31) , relevant for BBN, is expected to be as large 
as few percent. To reach accuracy of the order of one percent in rate evaluations it is 
therefore necessary to correct Born rates for both radiative and Coulomb effects. This last 
contribution, however, is only present when both electron and proton are present in the 
initial or final state, since it can be viewed as the electromagnetic rescattering of the two 
charged particles. In this way it is straightforward to show that no Coulomb corrections 
J'{E) are present for the channels e"*" + n ^ Vg + p. The results reported in ( |3.23| ), 
( |3.28|) can be applied to the general formula for the decay rates ( 3.10 ) by rewriting the 



dependence on the electron spectrum end-point Em in terms of neutrino energy. In fact. 



with the same notation of Figure 3.1 we have that E^ = Ei + E^. Then, in general the 
correction factors depend on both the electron energy E^ and neutrino energy Ei, which 
are both integration variables. 



Indicating with dTs the uncorrected decay rate in ( 3.1C| ), the corrected one is given 



■^For completeness we have included also the corrections coming from the finite nucleon mass changing 
the squared matrix element, which are considered in the next chapter. However, these corrections are very 
small. 
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Figure 3.4: Radiative corrections to the Born rates for the processes in ( |3.1| ), ( |3.2| ), (| 

by 

dVr ^ giEuE3)HE3)dTB , (3.32) 

where the factors are the same ones reported in ( |3.23|) , (|3.28| ) but with Em substituted by 
El + E3, and the factor T{E^) does not apply to reactions (c), (d). 

Equipped with the corrected formula in ( 3.32D , we may now numerically evaluate the 



Born rates for the reactions in (|3.1|), (|3.2D, (|3.3|) . In Figure 3.4 we show the relative 



difference between the Born rates and the radiatively corrected ones, namely the quantity 

AF. = ^^y^ . (3.33) 

Note that at sufficiently high temperatures the corrections for a given process and the 
inverse one coincide. Apart from the /3-decay reactions, whose corrections are practically 
constant with temperature, for the scattering reactions the considered radiative corrections 
are relevant especially at low temperatures. 

In our BBN calculations reported in chapter |6| we have also included an overall rescal- 
ing factor /t^^^ for the reaction rates allowing the theoretical prediction for the neutron 
lifetime to be fully compatible with the experimental value. This is a standard procedure, 
which allows to overcome the problem of a precise determination of the coupling constant 
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for weak interactions involved in BBN, expressing this overall factor in terms of the ex- 
perimental value of neutron lifetime. We stress, however, that it is worth rescaling the 
rate only after all known corrections have been included, since this increases the accuracy 
of the prediction. In this way the ansatz that the residual correction be an overall factor 
and not, as it is reasonable to expect, a function of leptons energies, may introduce errors 



less than 1%. Although not included in Figure liA (it would cause a trivial overall shift of 
the curves), we have also rescaled the corrected rates in ( 3.32| ) by the factor 

j-th 

1 + '^r = 4W - 1-008 , (3.34) 

Tn 

representing an energy independent constant correction which should be included to re- 
produce the experimental value for neutron lifetime at zero temperature and density. 
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Chapter 4 

Finite nucleon mass corrections 



In this chapter we calculate the corrections to the Born rates for the n ^ p reactions (|3.l| ) , 



( |3.2D , (3.3) induced by relaxing the assumption of infinite mass nucleons |6^. Such an 
assumption was made (implicitly or explicitly) in several points of the calculations done 
in the previous chapter, which will be now subject to critical examination. 

The main effect of considering finite values for the nucleon masses are due to a non 
vanishing nucleon velocity (and then recoil effect) and to the presence to effective extra 
nucleon interactions such as weak magnetism. 

Recoil effects have been neglected during the calculation of Born rates in the evaluation 
of the spin summed squared matrix element (Eq. ( p. 12 ) is valid only for fixed nucleons) and 
in the calculation of phase space integrals. The last point (kinematical changes) comes out 
from the fact that a finite nucleon velocity modifies the energy-momentum conservation 
relations and then the phase space for the given process changes. In particular, both the 
statistical distributions (Fermi functions) and the integration limits change with respect 
to the Born case (note also that Eq. ( |3.9D is valid only if one neglects recoil effect in 



Sspins l-^P)- The most practical implication of considering recoil effects is that it is no 
more possible to write the reaction rate as a one-dimensional integral, as in Eq. ( |3.16D . 

However, recoil effects are not the only ones involved in the finite nucleon mass cor- 
rections. Up to the same order of approximation we have also to consider modifications of 
the effective nucleon weak current due to weak magnetism, and, in general, to interaction 
terms due to the scalar and pseudoscalar couplings of the nucleons. 

In the following section we firstly calculate the dynamical changes in the squared 
matrix element induced by weak magnetism, scalar and pseudoscalar interactions. Then, 
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in the subsequent section, we will examine kinematical changes induced in the phase space 
structure of the reaction rates. At first order in T /Mj^ these two corrections can be treated 
independently. 



4.1 Corrections to the transition amplitude 



For definiteness let us consider the process shown in Figure ( |3.1| ). At first order in I/Mat, 
the transition amplitude can be written as 



M = —Up{pi)O^Un{p2)Ue{pz)l^{l-l5)Uu{Pl) , (4.1) 



where the effective nucleon-nucleon weak coupling is given in general by ||^ 

= if,{Cv - CaI^) + ^ Vi— o-^i' + h% + fpslbQf, ■ (4.2) 

Here denotes the momentum transfer to the final nucleon and /2, /a, fps are the 
anomalous weak charged-current magnetic moment, scalar and pseudoscalar couplings of 
the nucleon, respectively In general, the couplings Cy, Ca, /21 /si fps are form- 

factors, which all depend on q'^. However, at the relevant energy scales for the considered 
processes, this g'^-dependence can be neglected, being of higher order than l/M^. 

With some algebra, one can evaluate the squared modulus of ( [4.1[) summed over all 
spins. It is very useful, for references to the other processes, to write the resultant expres- 
sion in terms of the relativistic invariants s, 1 With the same notation of section ( |3.1| ) 
we obtain the result reported in appendix ^ The expression in ( PH) holds for all the SIX 
processes in (U), (U), (U) . 



Let us note that in the formula (EJ) we have retained only the leading terms, namely 
those corresponding to first order in the couplings /2, /s and fpg. Aa a further simplifi- 
cation we can drop out all the terms coming from scalar and pseudoscalar interactions, 
which are indeed very weak (in the following we then use /a = fps = 0) 



In Figure 4.1 we plot the "zero temperature" cumulative corrections, i.e. the cor- 
rections to the transition amplitude considered here plus the radiative QED corrections 
evaluated in the previous chapter; for the process ^ e~ +p we have AFr = F/j — F^ 



^The third invariant u has been eliminated through the relation u — — s —t + Mf + A/| + Mf + M| 
which follows from energy-momentum conservation. 
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Figure 4.1: "Zero temperature" radiative corrections to the Born rates for the processes 
in (|3.1|) , (|3.2| ), (^) (see text). We use the same notations of chapter ^. 
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ri?(fe + n ^ e +p) 



l + 6r 
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n ^0 



dEi j dEs \M\^{Ei,E3) 

^'^"f spins 



X g{Ei,E3) T{E^) F,{Ei) [1 - Fe{Es)] 



(4.3) 



where 



Einf 
Esup 



(M„ + ^i)(M2 - m2 + m2 + 2^1 M„) - 2^1^ 



2M„(M„ + 2^i) 
(M„ + ^i)(M2 - m2 + + 2EiMn) + 2Eii 



2Mn{Mn + 2Ei] 



M2 - M2 - m2 + 2£;iM„ 



4m2M2 



(4.4) 

(4.5) 
(4.6) 



The reason for this name is that in the hmit of zero temperature these corrections do 
not vanish, differently from what occurs to the other corrections considered in this thesis. 
Note that the constant shift correction ndr/a (see Eq. (|3.34D ) has been subtracted in 
order to show the pure radiative and Coulomb effects and finite mass corrections. As we 
will see, these corrections turn to be the most relevant ones for the Born rates. 
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4.2 Kinematical corrections 



Apart from the dynamical corrections to the squared matrix element calculated in the pre- 
vious section, there are also other corrections to the reaction rates which are of kinematical 
nature, arising properly from the finite nucleon velocity. While dynamical corrections are 
due to the fact that for a slightly moving nucleon the collision probability is higher than 
for a fixed one, kinematical corrections involve changes in the energy carried off by leptons 
due to nucleon motion. 

It is important to note that the reaction rates (per incident nucleon) intervening in the 
Boltzmann equations for the calculations of the elemental abundances have to be evaluated 
in the reference frame of the comoving volume (radiation rest frame) . 

For definiteness let us consider the reaction + n ^ e~ + p whose exact expression 
for the rate is given by 

T{Ue + n^e + p) = - TTT-To TTET ^n(^2) 



n J (27r)3 2E2 ^ ^ [J (27r)3 2E4 (27r)3 2Ei (27r)3 2E3 
{27r)U\pi+p2-p3-P4) E l^^l'f F,{Ei){l-Fe{Es)){4.7) 

spins J 

(as already discussed, we approximate 1 — Fp{E/^) ~ 1), where 

is the incident neutron number density. All the quantities in (4.7) are evaluated in the 



comoving volume reference frame; we denote with a prime the corresponding values in the 
neutron rest frame (in the infinite nucleon mass limit the two reference frames coincide). 
If the neutron motion lies along the z-axis, we have 

(:j-(s)^(;.7)(oO ■ 

where u = P2z/^h is the neutron velocity and 

-y = ^ ~ 1 + !!! = 1 + (4 10) 

= : ~ n = — — . (4.11) 

For the leptons we then obtain 



„2 
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where is the lepton momentum component along neutron motion. In general for the 
lepton energies we can write 

E,, E[, + ^^^^ + ^E[, , (4.13) 



Now, let us note that even if F in ( [4.7^ has to be evaluated in the comoving volume 
reference frame, the quantity in {...} is a Lorentz invariant which takes the same value in 
any reference frame. We choose to evaluate this quantity in the nucleon rest frame. In 
this way, the kinematical corrections involve only the quantities outside the braces 0. The 
correction to the phase space term is then easily calculated 

5(F,(1 - Fe)) = F,{E^) (1 - F(F3)) - F,iE[) (1 - F{E',)) 

^ F^{E[) (1 - F{E',)) i^^F{E',) - ^ (1 - F,{E[)) + 

1 /P2-P3V 



2 ^(^^) (1 - 2F(F^)) + 

H™1f(F^)(I-F.(F1)) + 

1 /P2-Pl\^ 



+ -^^E'.FiE'.) E\ il - FJE\))\ . (4.14) 

It is here appropriate to approximate the neutron distribution function with a Boltzmann 
function 

1 E2 Mo P2^ 

Fn{E2) = — ^ e'— ~ e-~ e ^"'2^ , (4.15) 

e~ + 1 

thus the correction to the rate for the considered reaction is 

1 f d^p2 1 '^■'2 H2L f f d^p4 1 d^pi 1 



6T {ue n —>■ e p) ~ — / - — — - ere 



n J (27r)3 2F2 [7 (27r)3 2F4 (27r)3 2Fi 

^ ^ i^^)' S\pi + P2 - P3 - P,) J2 HFu (l - F,)) . (4.16) 

y ' spins J 

Since the expression in {...} does not depend on neutron momentum (see the above dis- 
cussion), it is very convenient to perform first the integration over this quantity. First 

■^Note that even if the quantity in {...} is an invariant, nevertheless the integration limits, in general, 
change their values allowing changes in the phase space due to neutron motion. However, we only want to 
calculate the leading kinematical correction contribution, thus we have to consider only corrections coming 
from the terms outside braces, while the quantities inside them (including the integration limits) have to 
be calculated in the infinite nucleon mass approximation. 
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order terms in ( 4.14 ) do not contribute, since 



d^P2 1 



M2 







(27r)3 2M2 

for symmetry reasons. Instead for second order terms we have 



(4.17) 
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(4.18) 



hence Eq. ( [4.16| ) now becomes 



(5r(z/e n ^ e p) ~ y 



d3 



P4 



1 d^pi 1 (i3p3 1 



(27r)3 2^4 (27r)3 2^1 (27r)3 2^3 
(27r)^5^(pi+|52-P3-P4) E 1^1' 



(4.19) 
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(5$ 



T 
4M| 



1 P3 



^F(i?^) (1 - 2F(i?^)) - ^FiE',) (1 - F,(i?i)) + 



(4.20) 



and X^spins l-^P is given by the expression in ( p. 12 ). As in section (3.1), the integration 
over the proton 3-momentum can be ehminated by using the 6^{pi + P2 — Ps — P4), 
whereas the neutrino energy integration (for approximatively massless neutrinos) through 
the energy (5-function. Also integration over angles can be analytically performed (only 
the one over electron- neutrino angle is now non trivial). The result is 

T 



ATKi^e n ^ e p) 



gM4 + 3c' 
2^3 



2M„ 



J dp3 pI Q2 e{Q) F,{Q) (1 - F{E^)) ■ 



^(1 - F.(Q))(1 - 2F,{Q)) - ^F{E,){1 - 2F{E,))+ 
3^F(i?3)-3|-(l-F.(g))-3 + 



+ 



Cy + 2>C\ 



STT^Es 



F{E^){1 - F,{Q)) 



(4.21) 
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where we have used the same notation of section (3.1). The above expression for ATxi^'e + 
n — > e~ +p) cannot be extended to all reactions of (^Tll), ( |3.2D , ( |3.3[ ) by simply using the 
substitution rules of Table |3.1[ This is due to the different form of the statistical factors 
involved in the corresponding expressions of (4.7). However, it is easy to find that the 
factor in square brackets must be replaced as follows. First of all, note that for the reaction 
i^e +p ^ e'^ + n this factor is the same as in ( [4.21 ) . Differently, for + p ^ + n and 
+ n ^ Ve + p, one has 



|^F,(Q) [1 - 2F,(Q)] + [1 - F{Es)] [1 - 2F{E,)] + 3^F,(Q) 



-3 - 

-3^ [1 - FiEs)] + 



ri 2 /-» 2 




2P3^Q 

STT^E^ 



F,{Q) [l-FiEs)] 



(4.22) 



while for n ^ e + Vg + p 

-3 - ^F,{Q) [1 - 2F,m - ^F{E,) [1 - 2F{E,)] + 3^F,{Q) 



+ 3^F(£;3) 



Finally for e + Ue + P 



r< 2 r< 2 

:cv^+3c7 



2P3'Q \ 

3TT,E3 J 



F{Es)F,{Q) 



(4.23) 



n 



-3 + 1^ [1 - F,m [1 - 2F,{Q)] + 1^ [1 - F{Es)] [1 - 2F{E,)] " 3|- [1 - i^.(Q)] 



3^ [1 - F{E,)] 




n 2 n 2 




2P3^Q 

3TT^E3 



[l-F{Es)] [l-F,m 



(4.24) 



We report the results in Figure 4.2 in the form of the ratio AT k / {a / b for comparison 
with the other corrections which are of order a. Note that in the entire temperature 
range relevant for nucleosynthesis the dominant corrections are positive and apply to the 
scattering reactions (b), (d). 
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Figure 4.2: The ratio l^T k / {ot / Ty)T b ■ Dashed hnes correspond to negative values and 
labels refer to reactions in (^.iD , (|3.2D , (|3.3D . 
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Chapter 5 

QED thermal radiative corrections 



The reaction rates for the processes in (3J), (|3.2|) , ( p7^ ) strongly depend whether they 
take place in vacuum or in a heat bath at given temperature T. The most important 
contribution comes from the density of states in the integration over phase space. This 
has been already taken into account in the previous chapters for the calculation of the 
rates in the Born approximation and of the finite nucleon mass corrections. However, 
there are also a number of effects [^], [^2| - [^] induced by a T 7^ background which 
have to be considered; they are usually referred to as thermal radiative corrections. 

Since the temperature range of interest for BBN is around 1 MeV, the most significative 
thermal radiative corrections to the processes in (p.lj), ( |3.2D , ( |3.3| ) arise in the context of 
QED. First of all, if the particles entering in the given process propagate in a heat bath, 
they in general acquire an effective thermal mass induced by coherent interactions with 
the medium (in particular with the photons in the bath). The Feynman diagrams to be 



considered in this case are reported in Figure 5.1 for the process (^). 

Moreover, wavefunction renormalization caused by T 7^ background has to be taken 

into account as well. Finally, there are vertex corrections mediated by the exchange of a 



photon in the bath; they are sketched in Figure 5.2. These corrections contain infrared 
divergences, which have to be eliminated via the inclusion of corrections coming from the 
spontaneous and induced emission and absorption of thermal photons, as described by the 



diagrams in Figure 5.3. 

All these effects will be here considered at order a. While we have to evaluate the 



square of the diagrams in Figure 5.3, for the mass shift and vertex corrections only the 



interference of the diagrams in Figure 5.1 and 5.2 with the Born graph in Figure |3.1| will 
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give contribution of order a. Note that we have considered only electromagnetic interac- 
tions of electrons and protons interacting with the particles in the plasma through their 
electric charge; this gives indeed the main contribution to the effect we are considering. 
Moreover, since these effects already contribute as corrections to the Born rates, they will 
be calculated in the infinite nucleon mass approximation. In particular we will adopt the 
nucleon rest frame. 

Finally, the modified dispersion relation for and photons also call for corrections 
to the thermodynamic quantities (energy density, pressure and so on) ||8^ calculated in 
chapter ||, which are relevant for the BBN itself. 

All these effects will be considered in this chapter; in particular in the following four 
sections we focus on the thermal radiative corrections to the Born rates for the reactions in 
( |3.lD , (|3.2| ), , while in the last section we consider finite temperature QED corrections 
to the relevant thermodynamic quantities. Note, however, that since we are looking only 
for the main contributions coming from these radiative corrections, we can safely assume, 
throughout this chapter, a negligible electron chemical potential, which is usually, in fact, 
a small quantity. 

Our calculations will be carried out in the framework of the real time formalism (RTF) 



[ p6|| of the finite temperature and density quantum field theory [p7| , 8S], an account of 
which is given in appendix ^. 

5.1 Mass shift correction 

Due to interactions with the particles in the plasma, an electron (or a proton 0) acquires an 
effective mass depending on the temperature of the bath. Therefore, the mass parameter 
entering in the expressions for the decay rates should be substituted with its effective 
value, which has been explicitly calculated in the appendix ^. 

We focus on the reaction Vf. + n — > e~ + p (see Figure |5.1| ) , whose Born rate is given in 
Eq. ( p. 161 ); in this expression, according to ( |F.33| ) (see appendix^), we have to substitute 

^Since electrons and protons, at a first approximation, have equal electromagnetic interactions, the 
self-energy terms for these two particles are equal. However, the mass shift correction for protons is much 
less important than that for electrons because of their large (vacuum) mass. We are only interested in 
leading terms and so we disregard proton mass shift corrections. 
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^ + fi and, correspondingly, 

Q ^ Q + , (5.1) 

in all quantities appearing in §11). With the same notation of the appendix ^, the 
parameter // is given by 



^l^H + l. r dk-F(E) + r dk-F(E) (log A -log 5)1 

vr \ 3 ^3 ^3 Vo E ^ ' ^ 2Esps Jo E ^ ' ^ ^ ^ 'j 

(5.2) 



with E = \Jk^ ^ m^, E-^ = + and A,B are reported in (|F.26 ). The mass shift 
correction we then obtain is 

AFm {i^e +n^e- +p) = ^^^3 ^ j dp3ple{Q) ■ 

■ {{Q + fif F,{Q + f,) {l-F{Es + fi)) - Q'F.iQ) (1 - F{Es))} , (5.3) 

Expanding in the variable fi, at first order in a we finally get 



ATMii^en^e-p) = ^ ^ • J dp^plQ" ^i9{Q) F,{Q) (1 - ^(^3)) • 

2 F(^ _ 

Q T T, J ^ ' 



Similar expressions for the other reactions in {3A), (^), ( |3.3D can be easily found by 



using the substitutions of Table |3^ and replacing the factor 



' + ^ - ^ ) , (5.5) 



in the above expression as follows. This factor applies to the reaction + p ^ e"*" + n as 
well. For e~ + p ^ + n and e"*" + n — > I/g + p, we have instead 

2 l-EiE,)^E^^ ^ ^^^^^ 



(5.7) 



Q T n 
while for n — > e~ + I7g + p 

2 F{E^) F,{Q) 

Q T 



and e + Ue + p ^ n 



2 l-F{Es) l-F^iQY 

Q ^+ n ' • ^^-^^ 
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Figure 5.4: The ratio l!^T m / {oi / ti)T b ■ Dashed hnes correspond to negative values and 



labels refer to reactions in (p.lD , (3.2), (3.3) . 



In Figure 5.4 the ratio l^T m / {ol / 'k)T b is reported as a function of the (photon) temperature 



T for the six reactions in (3.1), ( |3.2D , ( p.3| ) . Note that mass shift corrections are negative 
for all these reactions and decreasing (in absolute value) with decreasing T for (a), (b), 
(e), (f) while they are almost constant for (c), (d) in the temperature range relevant for 
nucleosynthesis. 

5.2 Wavefunction renormalization correction 



The same diagrams in Figure |5.1| leading to mass shift corrections also contribute to wave- 
function renormalization corrections. For definiteness we consider in detail the reaction 
( |3.lD . As studied in appendix ^, the amplitude is again given by (|3Tll ), but when we 
compute X^spins l-Wp we have to replace the projectors with the renormalized ones, ac- 
cording to ( W^ - With the same notation of chapter |3| we obtain the following expression 



for the correction to the decay rate: 



/Wwii^e +n ^ e- +p) = ^^^^ J dp3plQ^e{Q)XF,{Q)il-F{ps)) (5 
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where the parameter A is reported in the appendix ^. Note that, in the non-relativistic 
hmit for the nucleons, which we here adopt, the second term in ( |F.42 ) contribute to the 



correction AF^i/ with a term which vanishes after the angular integration is performed. 
From relations it is easy to show that 

A = ;^ {- {f^ + ) + (ff + f/') + me (c^ + c'^)} (5.10) 

The quantities f^, ff can be immediately read off by Eqs. ( [F.15D , ( |F.lg| ) by substituting 
Pq = £'3. Furthermore T^', Tp', c'^, c'p may be calculated by differentiating Eqs. ( |F.14| ), 
( IF.17D , ( |F.16| ), (|F.19| ) with respect to 

Po prior the integration in dx is performed. After 

some algebra, the final result is 

TT J k 6Esp3 Es - ps 

- ^ Z fdk^ F{E) iiE3 + E)lnA - (E^ - E) In B) + 
2tt Esp^ J E 

2a /■ „ 1 ^ ^ , , 



Note that the first term in (5.11) is infrared divergent. This term, together with another 
similar divergent one coming from considering vertex corrections (see the following sec- 
tion), is cancelled by an opposite term provided by the photon and emission rate discussed 



in section 5.4. The collinear divergence still present in ( [F.38D for p^ = k, together with 



a similar contribution from the vertex correction is also compensated by bremsstrahlung 
diagrams. 

The extension of ( ^ ) to the other processes in (0), ( |3.2| ), (|3.3|) can be straightforwardly 
obtained by using Table |3.1| . 

The results for ATw/{a/'n')TB are plotted in Figure ^.5| . Also these corrections are nega- 
tive for all the six reactions in (3.1), (p.2|), (3.3) , and are monotonically decreasing with 



decreasing T. 

5.3 Vertex correction 

The vertex correction to the Born rate for the process + n — > e~ + p is given by the 



diagram in Figure |5.2| , whose amplitude is 



Mv = y" ^^^Mp(?'4)(ie7^)5'p(p4 + A:)7^(Cv-C^75)nn(p2) X 

X iDp^{k)ue{p3){ie-f'')Se{p3-k)'y^{l--f5)u^^{pi) , (5.12) 
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Figure 5.5: The ratio ISTyi/ / {a/-K)T b- Dashed hnes correspond to negative values and 



labels refer to reactions in (p.lD , (3.2), (3.3) . 



where Sp (F = e,p) and Dp^^ are the real time formalism fermion and photon propagators 
respectively, reported in the Appendix^. After some algebra, and using the Dirac equation 
for the proton. My can be cast in the form 



where 



(5.13) 



= He 7^(1- 75) , 



and 



2 I ^k + p3)mp4+P3 + ky-M^) 
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S{k'^)k^k'' Adik"^ -ml)F{ko){p3 + k)^ip3 + kY 



P4,-kp3-k {k + P3y {{Pi + P3 + k)'^ ~ M"^) J 

The order a vertex correction to the Born rate is provided by the interference term of the 



diagram in Figure 5.2 with the tree level diagram in Figure |3.]j , 

MvM^ + M^Mb , (5.14) 



where Mb is the Born expression (3.11) 



Gf 
V2 



Mb = i^Af.B'' . (5.15) 



After some lengthy calculations, the final result for the vertex correction we obtain is 

ATviiye + n^e-+p) = ^— ^ / dp3 P3 dk k 0{Q) 

ivr-^ vr Jo Jo 

.FM)l^-Fmi^{^^o,A^^<.,B--^}. . (5.16) 

Note that we have only reported the non infrared divergent part (as for ATyy, the infrared 
one comes from the Bose terms only, and will be cancelled by the bremsstrahlung contri- 
bution) 0. 

The expression in ( 5.16| ) is extended to the other processes in (|3.1|), (p. 2]), ( |3.3| ) by using 



the Table 3.1, and the plot for these corrections is reported in Figure |5.6| . Note that for 
the scattering reactions (a), (b), (c), (d) the vertex corrections are positive in the relevant 
temperature range, while for (e), (f) the dominant contribution is negative. 

5.4 Photon emission and absorption 

In the primordial plasma, in which an electromagnetic component is present, spontaneous 
and induced photon emission as well as photon absorption take place, so that when con- 



sidering a given process such as the one corresponding to the diagram in Figure 3.1, the 



bremsstrahlung processes depicted by the diagrams in Figures 5.3 must be included as 



well. These also cancel the infrared divergences due to the radiative diagrams of Figure 



^Note that the result for the vertex correction quoted in [[72| as Eq. (11) has a missing factor 1/E. 
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Figure 5.6: The ratio ISX^y jip.j'K^V b- Dashed lines correspond to negative values and 



labels refer to reactions in (p.lD , (3.2), (3.3) . 



|5.2| corresponding to wavefunction and vertex renormalizations. 
The matrix element for the emission processes in Figures is the following: 



V2 



{Up{pA) l^{Cv - CaIs) Un{p2) Ueips) (^67") Se{p3 + k) 7^(1 - 75) U^^{pi) + 



Up{pA) (ie7") Sp{pi + k) Y{Cv - C475) Un{p2) Ueips) 7^(1 - 75) u^M} 4{k) 



while a similar expression, M^, holds for the absorption processes. In order to simplify 
the compuatation of the spin summed squared matrix element, one can subtract from 
the beginning the infrared divergent terms coming from the propagators. They are all 
proportional to the Born X^spins ^-nd will cancel the divergences in AFyi/ and AFy. 
Note also that, for kinematical reasons, only the T = part in the propagators contribute 
to the amplitude. 

After some extensive calculations, we obtain the following finite contribution to the rate 



AF., = AF^ + AF^ 



AF^(fe + n^e + p) 



X [l-F{Es)] 



2 

A u a 



27r3 



dp3 rdk4-B{k) 

vr Jo Jo -C/3 



2^, 



J^log (^l±Pl 
kp3 \E3-p3 



Ql + Ql- 2Q2 
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Figure 5.7: The ratio l!^T / [a / ii)T b ■ Dashed lines correspond to negative values and 
labels refer to reactions in (p.lD , (3.2), (3.3) . 



2i-^ log -^-^ 



Ql +— log 



E3+P3 



E:t 



P3 



(5.17) 



where Ql = {Q±kfF^{Q±k) 9{Q±k), = Q^F^{Q) e{Q) and 5 = 1 Q. 

For the other processes in ( |3.lD , (^) , (|3.3| ) we have to use the Table |3]l| and set ^ = 1 for 

Ve + n <-> e~ + p and n <-> e~ + I'e + and ^ = for e"*" + n <-> I^g + p. 



The results for H^T ^ / [a / ti)T b are reported in Figure 5.7. It is easy to realize, from 
this figure, that photon emission and absorption contribution to the thermal radiative 
corrections for neutron decay and the inverse process is the dominant one (and positive) 
in the temperature range relevant for nucleosynthesis. This is simply understood since the 
inclusion of the process 7 + n — > e~ + Vf. + p greatly increases the neutron decay rates, 
otherwise strongly suppressed by phase space. 



■^Also here we must note that our result corrects Eq. (13) in Ref. |7^] for the photon emission and 
absorption correction previously calculated. 
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5.5 Corrections to the equation of state 



The temperature dependent effective mass of electrons and positrons in the primordial 
plasma affects the calculations of the thermodynamic quantities like particle number den- 
sity, energy density, pressure and entropy, which directly enter in the evolution equations 



for the light element abundances (see section 2A), as well as in the weak decay rates as 



discussed in section 5A. In addition, also the induced effective photon mass has to be 
taken into account. The modified dispersion relations of e^, 7 have been calculated in 
appendix ^ and they must be used (instead of the vacuum dispersion relations) in the 
evaluations of the quantities in (|L4^-(|L5^). 

The relevant quantities to be considered are the e^, 7 energy densities, entering in Eqs. 



( CT) and ( ^) , and pressure, which enter directly in the Eq. ( p.6[)|) for entropy density. 



Corrections to electron and positron number densities are also present in general, but, 
since in the electric charge conservation equation ( p. 461 ) only their difference is relevant, 
this results into a correction to the electron chemical potential. This is, generally, a small 
quantity and the corrections to it can be usually neglected. 

The order a corrections to the photon energy density and pressure are easily calculated 
by inserting Eq. ( [F.57 ) in the Eqs. ( 1.49|) , ( 1.50|) specialized to photons; with simple 



manipulations we obtain 

dp - 

/o 



Spy ^ T'- / dp^-F,{E) , (5.18) 

IT Jo E 



Sp-r - -n^P-y ■ (5.19) 



The expressions for order a corrections to energy density and pressure are deduced by 
substituting Eq. ( |F.33D into ( 1.491 ) , ( [L.50D specialized to ; we found (pe = Pe- + Pe+ > 

Pe = Pe- +Pe+) 

5pe - ^ dpp" p[l - -Fe{-E)^ F,{E) , (5.20) 
Spe ^ ^ dp^p[l + -F,i-E)jF,{E) (5.21) 



(E = ^p'+mi). 

The corrections to the energy density and pressure of the electromagnetic plasma {5p = 
{m/-iT^T'^){5r^ + 6pe), (90/7r2r4)(,^p = 5pj + Spe)) are plotted in Figure |5]|. As we can 

75 




Figure 5.8: Radiative corrections to the energy density (solid line) and pressure (dashed 
line) of the electromagnetic component of the primordial plasma. 

see, its contribution is of the same order of magnitude of the other thermal radiative 
corrections. 

As briefly mentioned above, corrections to the energy densities directly influence the 
expansion rate of the Universe (see Eq. ( |2.44| )) and eventually the primordial abundances 
through a change in the freeze out temperature. Moreover corrections to p and p also 
translate into corrections to the entropy density according to relation ( |1.79| ). On one 
side, this traduces into a modification of the time-temperature relation or, to be more 
specific, of the evolution equation for T ( p.60| ). On the other side, the entropy of the 
plasma is transferred to photons when pairs disappear, while neutrinos do not benefit 
of this entropy releasing since at this epoch they are decoupled. Hence, this changes the 
neutrino-to-photon temperature ratio, which directly affects the weak rates. We have also 
calculated this effect, but we found that the modification oiT^/T is quite negligible in the 
range relevant for BBN, much smaller than the precision of our goal. 
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Chapter 6 

Calculations of Big Bang 
Nucleosynthesis. Results 

In this thesis we have studied in details the fundamental processes underlying primordial 
nucleosynthesis and performed a complete analysis of the physical approximations that 
may be used as well as of corrections which should be included to reach the 1% reliability of 
the estimates. The relevance of this study and, in particular, of the precise determination of 
the primordial light element abundances, has been already pointed out in the Introduction. 
In this final chapter we summarize the main results obtained, give the accurate prediction 
for ^ife and outline the basic ingredients of forthcoming numerical code computing light 
element abundances, which is currently in progress. 

6.1 Primordial nucleosynthesis: cosmology and particle physics 
inputs 

As discussed in chapter ^, big bang nucleosynthesis is the result of two competing factors: 
nuclear and subnuclear processes between the particles in the primordial plasma, whose 
effect is to keep the species in thermal equilibrium, and the expansion of the Universe, 
which tends to hinder this equilibrium. It is then clear that the physics involved runs over 
cosmology (and, in particular, thermodynamics in the expanding Universe) and particle 
physics. Here we mainly summarize the basic inputs for studying primordial nucleosyn- 
thesis as well as the physical assumptions and approximations we used. 

Standard cosmology is based on the observed homogeneity and isotropy of the Uni- 
verse. These directly lead to the cosmological equations ( |1.28| ), ( |1.29| ) in the approximation 
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of the Universe as a perfect fluid, to which must be added the equation of state ( 1.30 ) for 
the matter fields present in the Universe must be added. Primordial nucleosynthesis takes 
place when the curvature term k/E? in ( 1.2^ ) is completely negligible, so that it is usually 
dropped out from the subsequent considerations. 

From the cosmological equations, the energy conservation law ( p.. 27 ) follows, as well as 
entropy conservation ( |l.78| ) for species in kinetic equilibrium. The total energy density 



and pressure, for species in equilibrium, are given by ( 1.84 ), ( 1.85 ). From these it follows 
that non relativistic matter (primarily baryons) gives a small contribution to the evolu- 
tion in the RD era (in which BBN takes place), since its energy density is exponentially 
suppressed. This approximation will be used in several occasions. 

The approach and the departure from thermal equilibrium is fully described by the Boltz- 
mann transport equations, which involve the interaction rates of the particles present in 
the plasma. It is found that when a ultrarelativistic or non relativistic specie is completely 
decoupled from the heat bath, its distribution has the same form as that of an equilibrium 
one, but with a temperature parameter scaling, as the Universe expands, according to 
Eqs. ( 1.96| ) or ( |1.99| ). For such decoupled species, entropy is then separately conserved. 
Further ingredients come from particle physics; let us focus on the nucleosynthesis era, 
when the relevant particles present in the plasma are photons, neutrinos, electrons and 
positrons, baryons. Photons and pairs form the background electromagnetic plasma, 
and their equilibrium is established by electromagnetic interactions. Baryons are main- 
tained in kinetic equilibrium by electromagnetic and strong nuclear interactions, while 
neutrinos by weak interaction processes only. Chemical equilibrium in the electromag- 
netic component of the plasma is established by the annihilation reactions e^e~ ^77, 



while the weak reactions in (3^), (3^), ( p. 31) are responsible for the nuclear chemical 
equilibrium. Departure from chemical equilibrium happens as follows. The electron neu- 
trino decoupling temperature is calculated to be T^^ ~ 2.3 MeV, while for v^^fr we 
have '^^ ~ 3.5 AleV . In nucleosynthesis calculations it is assumed that neutrinos are 
completely decoupled from the background plasma before e^e~ annihilations. Neutrino 
decoupling, however, occurs very close to e^e~ annihilation, so that some residual interac- 
tion with the thermal plasma can cause the neutrinos to be slightly heated by the resultant 
entropy release. Careful studies have been conducted on this subject [^, Ssl, demonstrat- 
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ing that neutrino decoupling is not an instantaneous process. Hence, the spectrum of the 
decoupled neutrinos deviates slightly from the Fermi-Dirac form, causing the effective neu- 
trino temperature to increase with momentum. However this increase accounts for only 0.7 
% (even at relatively high momenta), justifying the usual approximation of instantaneous 
decoupling. Furthermore, the electron chemical potential, which is constrained by electric 
charge conservation, as well as the baryon chemical potential, is usually supposed to be a 
small quantity. This assumption is practically justified a posteriori. Instead, the neutrino 
chemical potentials are unconstrained, but they are set to zero in the standard nucleosyn- 
thesis scenario. Note that a non zero chemical potential for electron neutrinos can alter 
neutron-proton equilibrium, as well as increase the expansion rate of the Universe, while 
chemical potentials for other neutrino types can only speed up the expansion. In general, 
the lowering of the n/p ratio at freeze out may be compensated for by the net speed up 
of the expansion rate, although exotic models exist in which large lepton numbers can be 
generated (see ||7) and references therein) Q 

There are no further assumptions and approximations regarding thermodynamics. 

Electromagnetic and weak interaction reactions between elementary particles are cal- 
culated in the framework of the electroweak Standard Model. In particular, the crucial (for 
nucleosynthesis) weak processes in (|3.1| ), (^), (^) have been here evaluated as follows. 
Born approximation is used as the reference for the reaction rates: it is based on tree level 
calculations of the probability amplitude in the infinite nucleon mass limit. We have then 
relaxed these assumptions. Firstly we reported the relevant zero temperature QED radia- 
tive and Coulomb corrections to the weak rates and then considered finite nucleon mass 
effects, namely weak magnetism, phase space modification and kinematical corrections due 
to the thermal motion of the initial nucleon in the comoving reference frame. The results 
for these corrections are plotted in Figures |4.1| , |4.2| . Furthermore, we have also studied in 
details thermal radiative corrections, explicitly induced by the finite density of the plasma. 
In particular, we have calculated electron effective mass and wavefunction renormalization 
effects, vertex corrections and thermal photon emission and absorption processes, and the 

^ These models are, however, constrained by BBN itself, since the introduction of non vanishing chemical 
potentials for neutrinos can substantially modify abundance predictions, destroying the agreement between 
theory and observations. In practice, allowed neutrino degeneracy cannot significantly altet the standard 
BBN model. 
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obtained results are reported in Figures 5.4, 5.5, j.t, 5.7. With reference to the above 



mentioned Figures, the results for each reaction channel can be summarized as follows. 

(a) Ue + n ^ e~ + p 

For the crucial BBN temperature range, 0.1 MeV < T < 3.5 MeV, the two main 
corrections to the Born rate come from zero temperature radiative and kinetic terms. 
The contribution AF/j is weakly depending on T and represents the dominant term, 
though for large temperature the kinetic correction AF;^ starts contributing signif- 
icantly. The two combined contributions correct the Born rate for a factor 6 9%, 
whilst thermal radiative ones are of the order of 1%. 

(b) e~ + p ^ i/e +n 

For this channel radiative corrections are dominant at low temperature, while the 
kinetic ones give a quite relevant effect in whole BBN temperature range and correct 
the Born rate for a factor varying from 1% at low T up to 3% for T = 3 4 MeV. 
The radiative contribution is quite rapidly decreasing with temperature, reaching 
large negative values. This leads to a partial cancellation between AT ft and AF^^-. 
Thermal corrections are dominated by bremsstrahlung contribution AF^ and can be 
as large as 2% of F^ for large temperature. 

(c) e"*" + n — > I7g + p 

For this process the radiative corrections have a behaviour quite similar to channel 
(b), though they are even more rapidly decreasing with temperature. This again leads 
to a partial cancellation between AFr and the positive monotonically increasing 
ATk- Thermal corrections are again mainly provided by photon emission/absorption 
and monotonically increase with temperature up to a factor 2% of the corresponding 
Born rate. 

(d) z7e + p ^ e'^ + n 

The radiative and kinetic corrections sum up to a factor of about 5% of the Born 
rate. Actually the opposite behaviour of ATr and AF^^ conspires to give an almost 
constant total correction in the whole interesting temperature range. Thermal effects 
are quite small, contributing for less than 0.5%. 
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(e) n ^ e + z^e + P 

For the neutron decay the radiative corrections are practically constant and give 
the leading effect for small temperature, T < nif.. For larger T the thermal photon 
emission/absorption processes rapidly become dominating over all other correction 
terms with a relative ratio to the Born rate as large as 10^. However, this large 
correction is weakly contributing to total n ^ p rate, since in the temperature range 
T > TTie the scattering processes {a)-{d) are largely dominant over decay and inverse 
decay. 

(f ) e~ + Ue +p ^ n 

Same considerations of the direct process (e) hold for inverse neutron decay. As 
for the direct channel, the kinetic contribution is negative down to temperatures of 
the order T ~ 0.2 MeV. The contribution of jAPi^l is however negligible for both 
processes, smaller than 1%. 



All results are summarized in Figure p.l[ where we have shown the total relative correc- 
tions in percent. For the processes (a) and (d) the correction is almost constant over the 
entire considered range for T, and of the order of 6 10% and 5 -r- 6%, respectively. The 
positive kinetic contribution soften the deep decreasing of the radiative terms for channels 
(6) and (c). Finally the large effect of thermal bremsstrahlung for neutron decay and 
inverse process (e) and (/) is particularly evident. The cumulative "zero temperature", 
kinematical and thermal radiative corrections to the n ^ p and p ^ n rates are instead 



reported in Figures 6.3, 6.4. We have shown the total rates T{n p) and T{p n) 
in Figure 6.5, while the total relative correction AF/r^ = {T — Tb)/^b^ in percent, are 
plotted in Figure |6.6| . The whole correction AF for n <^ p turn out to be a positive 
decreasing function over the whole temperature range relevant for BBN. The main con- 
tribution at low temperature for both total rates comes from the radiative corrections, 
while for T > 2 3 MeV kinetic contribution starts dominating. This is particularly 



evident by looking at Figures 6.2 and 6.3. While for radiative corrections the effect on 
n — > p total rate is larger than on the p — > n one, AF^^ shows an opposite behaviour. The 
competition of these two corrections is then responsible for the presence of the inversion 



points in Figure 6.6 at T ~ 0.15 MeV and T ~ 2 MeV. Finally, the order of magnitude 
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T/m. 



Figure 6.1: The total relative correction to the Born rates, in percent, are shown for the 
six processes in (|3.1|), (|3.2|), ( |3.3| ) . 



3 
<1 




1 



T/m, 



Figure 6.2: The zero temperature radiative corrections, for n — > p (sohd hne) and p ^ n 
(dashed line) processes. Finite mass contributions coming from phase space integration 
and weak magnetism are also included. The result is expressed in percent. 
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Figure 6.3: The kincmatical contributions to the n ^ p (solid hne) and p —> n (dashed 
line) processes, normalized to the Born rates, in percent. 




Figure 6.4: The whole thermal radiative corrections to the n — ^ p (solid line) and p ^ n 
(dashed line) processes, expressed in percent. 
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T/r 



Figure 6.5: The total rates T{n p) (solid line) and r(p — >■ n) (dashed line), including 
all radiative, finite mass and thermal corrections. 

of the pure thermal radiative corrections 

AFt = ATm + ATw + AFv + AT^ , (6.1) 

is sensibly smaller, but nevertheless they may contribute for a factor 0.2 0.4% at the 
freeze-out temperature T ~ 1 MeV. We have performed a fit of the numerical results for 
r(n — p) and r{p —>■ n). The fitting expressions are the following 

1 ^ / T \' 

1 ( qm \ ( T \} 

r(p-n) = ^exp^-^j , (6.3) 

where for n — >^ p 

ao = 1 ; ai = -0.0135609 ; a-i = 1.227825 ; 03 = -19.344104 ; 

04 = 85.9281797 ; as = -11.221606 ; ae = 14.4804529 ; ay = -3.082314 ; 
08 = 0.357466 ; og = -0.0200726 ; oio = 0.10049-10-3 ; an = 0.488642-10-^ ; 

ai2 = -0.237462-10-5 , ai3 = 0.359274-10-^ (6.4) 
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T/m. 



Figure 6.6: The total relative corrections AT/Tb for the n ^ p (solid line) and p ^ n 
(dashed line) processes, expressed in percent. 

while for p ^ n 

bi = 20.884200 ; 62 = -71.802238 ; 63 = 118.853956 ; 64 = -11.564594 ; 
65 = 45.057913 ; 66 = -3.316009 ; 67 = 0.2744236 ; 63 = -0.75654-10-2 ; 
69 = -0.37916-10-3 ; 610 = 0.201124-10-^ ; 611 = 0.791665- 10"^ ; 

612 = -0.640452-10-'^ ; 613 = 0.104891- IQ-^ ; g = 2.90377 . (6.5) 

The fit has been obtained requiring that the fitting functions differ by less than 0.1% from 
the numerical values. 

Finally, note that mass shift radiative corrections affect not only the weak interaction 
reaction rates but also the energy density and pressure of and photons in the Universe. 
We have considered the corrections to these quantities as well, and the results are reported 
in Figure 

Last but not least, we have to consider the strong nuclear processes intervening in pri- 
mordial nucleosynthesis. Unfortunately, there is no very reliable theoretical computation 
for the interested nuclear reaction rates, so the strategy adopted is quite different: exper- 
imental data on the processes involved are used. However, a certain degree of uncertainty 
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is introduced in this way, mainly due to the fact that the typical energies of the reactions 
taking place during BBN are quite different (too much low) from the available experimen- 
tal ones, so that it is necessary to take some extrapolations. Moreover, a suitable thermal 
averaging of the measured cross sections (usually made with a Boltzmann distribution) is 
required: 

<av>= / dE a{E) E e"^/^ (6.6) 

(M is the reduced mass). All this conspires to produce uncertainties in the obtained reaction 
rates as large as 10%. Fortunately, the impact of these on the predictions for light element 
abundances is extremely low (9^. This happens because, as we have already stressed, the 
crucial reactions for nucleosynthesis are those fixing the n/p ratio at freeze out (especially 
for '^He abundance), and such reactions are very well theoretically evaluated, as discussed 
in this thesis. All the sources of uncertainties in primordial nucleosynthesis predictions 
are, however, currently well controlled by both MonteCarlo [^T| simulations and semi- 



analytical analysis |9C] 



6.2 Prediction on the '^He abundance 

The conditions for BBN to take place were reached when the temperature in the Universe 
was in the range 10 0.1 MeV. Primordial nucleosynthesis is, in fact, the final outcome 
of the decoupling of the weak interactions which keep neutron and proton in thermal 
equilibrium and, shortly after, of the onset of nuclear reactions which start producing 
light nuclear species. At temperature T >> 1 MeV the neutron to proton ratio is order 
unity and starts decreasing exponentially when the temperature reaches the value of their 
mass difference. As the temperature decreases, however, weak interactions are no longer 
fast enough to maintain equilibrium and a substantial final neutron fraction survives down 
to the phase of nucleosynthesis. All neutrons become practically bound in ^He nuclei, due 
to the high binding energy per nucleon, which therefore represent the dominant products 
of BBN. While the '^He mass fraction is weakly depending on the baryon to photon 
ratio rj, it is strongly dependent on the neutron fraction at the nucleon freeze-out. An 
accurate theoretical prediction for the helium abundance, as well as for the other light 
nuclei produced during BBN, will be the subject of the next section; it can be obtained 
by using the standard BBN numerical code suitably modified to take into account all 
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the corrections considered in this thesis. Here we report on the results for the expected 
variation of the surviving neutron fraction Xn induced by the whole effects Ar(n ^ p). 
This allows for a simple estimate of the corresponding variation of Xijj^; the baseline 
value for this quantity obtained with the standard BBN code without the inclusion of any 



correction is |12] 



X^He = 0.2411 . (6.7) 
The neutron fraction is determined by the differential equation 

^ = ^[r(p-n)(l-X„)-r(n-p)X„] . (6.8) 

Writing X^ = X^ + where 5Xn is the correction induced by AT, we have at first 
order 

dX^ dt 



n 



J'Bip - n)(l - X") - TB{n ^ p)XOj , (6.9) 

^SXn = -^[{TBip^n) + TB{n^p))6Xn-AT{p^n) 

+ (Ar(p^n) + Ar(n^p))X°] . (6.10) 

Notice that the zeroth order abundance X^ has been defined as the one obtained by the 
Born amplitudes rescaled by the constant factor 961/886.7, which provides at tree level 
the correct prediction for neutron lifetime (see our discussion in section 3). The Born rates 
^b{p ^ "-) in ( |6.9| ),( |6!ToD are therefore rescaled by the same factor. Equations (|O|),( |O0D 
have been numerically solved using our fitting function for T{n <-> p) and a similar one for 
r_B(n <-> p), which we do not report for brevity. We found for the asymptotic abundance 
6Xn ~ 0.0024, with a relative change, in percent, SXn/X'^ = 1.6%. 

These results allows for a simple estimation of corrections to '^He mass fraction by means 
of the formula ( pO| ) 

X,H, ~ 2Xnitns) . (6.11) 

Note that X„(t ns) is evaluated at the time when nucleosynthesis begins, tns — 180 s; 
actually, by this time the neutron abundance surviving at freeze out has been depleted by 
/3-decay to 

Xnitns) ^ X^e-*-/^" , (6.12) 
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where r„ is the neutron hfetime. Hence the corrections to ^He mass fraction can be 
evaluated from 

5X4He ^ 2 5X„e-*""/^" . (6.13) 
and using the results for dX^ we find 

5XiHe - 0.004 ~ 1.6% . (6.14) 

The total correction, as mentioned in the previous section, is largely dominated by "zero 
temperature" radiative and finite nucleon mass corrections, which give a positive contri- 
bution to Xifjg^. 

6.3 Building the BBN code 

In this section we finally describe the way the primordial abundances of light nuclei can 



be numerically evaluated and discuss how to build of a BBN code - [62]. Such a code 
follows the evolution of these abundances, and primarily tracks p, n, D, T, ^He, '^He, ^Li, 
Li and ^ Be. The initial temperature can be chosen to be, for example, T = WMeV, 
and hence for the initial abundances we can take the nuclear statistical equilibrium values 

...(i|)-f.i(Jl)'-',-.|xr^e* (0.., 

(hereafter we use the same notation of chapter ^). In fact, at temperatures ^ IMeV, 
the nuclear rates are sufficiently high to keep all abundances at their equilibrium values. 
It is usually made the reasonable assumption that the elements are always maintained in 
kinetic equilibrium by electromagnetic and strong interactions, while this may be not true 
for chemical equilibrium; then the coefficients in the nuclear rates depend only on r] and 
T. 

The set of evolution equations relevant for nucleosynthesis, which is the main body of the 
BBN code, has been deduced in chapter |2[ Here, we report these equations for complete- 
ness: 



R dt 



= -3 4^ (6.17) 

UB dt R dt ^ ' 

88 



TT' UB qB 



dT 

n 



2 T^f{z) 
1 dR 

^l^~ir (Pl + P7 + Pe + Pe + QiT - Td) {pu + Pu) + Pb) 

R dt 



+ 



dpe I ^^d4>edYj_ 1 dR ^ 



dt 




+ ubY. (am, + -T 

dp. 



j,k,l 



NilNkl 



ij—>kl 



dT 

yN, yN, 
i j 



(6.18) 

dY^ 
dt 

(6.19) 
(6.20) 



As already mentioned, this set of equations does not include the evolution equations for 
neutrinos. Here, however, we have explicited the effect of not decoupled neutrinos on the 
involved quantities by means of the theta-function Q{T — Tq), Td being the decoupling 
temperature. The relevant thermodynamic quantities can be read off by Eqs. (2.51)- 
( |2.57 ), while qs is the baryon electric charge, as defined by equation ( 2.47 ). Note that it 
has been assumed that the electron chemical potential parameter (pe is small, so that all 
terms depending on (pg can be expanded at first order in this parameter. 
The other inputs are contained in the equations ( 6.2C| ), and are the weak interaction and 
strong nuclear reaction rates. The corrected weak n <-> p rates are taken from our fits 
in Eqs. (|6.2|)-( [6^ ). The complete nuclear reaction network is composed of 88 reactions, 
which we report in appendix ^ however, in many cases, one can also consider a limited 
number of reactions, introducing a negligible error in the computations. For brevity, we 
do not report the expressions for the nuclear reaction rates which can be found in ||& 

6.3.1 Numerical tricks 



The basics equations to be followed for evaluating the primordial abundances are ( 6.17] ) 
and ( |6.20| ). Eq. ( |6.16| ) can be directly substituted into ( |6.17 ), while the expression ( 6.18| ) 
for the electron chemical potential has to be used for calculating Pf>, Pe but also p^, p^ 
and p^ (see the previous footnote). Instead, by means of Eq. ( 6.19| ) we can translate the 



■^Note that, in general, both pe, Pe and p^, depend on the electron chemical potential, the last two 
quantitiex through the effective finite photon mass induced by the e^e~ bath. Furthermore, observe that 
from the entr opy coi iserv ation also p^, wo uld depend on </>£. However, this dependence is effective 
only in Eqs. (6.f6), (3.17), while in Eq. (3.19) the presence of the theta-function makes it non operative 
(neutrino entropy is related to the electron one only when they are coupled to the plasma. 
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time evolution into the temperature evolution of the interested quantities. 

For numerical reasons, it is better to turn the variable ub into the dimensionless 
quantity 

h , (6.21) 



which evolves more slowly with T than ns- Furthermore, instead of T, we use the variable 

(6.22) 

In terms of these news variables Eq. ( |6.17| ) and Eqs. ( 6.20| ) take the form 



me 
T 



dh 

dz 



(3p^ -p^ + 3pe -Pe + Z^ + Z§^ + fh Yj] 


H + hj:j[zAMj + f^ 




3 


Py + + Pe + Pe + ^Q{zd - z)pu + h J^j Yj ^ 







3/i 

z 



(6.23) 



dz 



4py + 4pe + i<d{zD - z)pu - z 



' dz 



3 [Py + p^ + Pe + Pe + ^@{zd - z)pu + hj^j Yj) H + hj^j [zAMj + 1) r 



z 



(6.24) 



with H denoting the dimensionless Hubble parameter H = H/rrif, 



H 



1 /2 

V^Si Py + Pe + P^ + h(^zMu + J2(^zAMi + ^^Y^ , (6.25) 

and ZD = me/To. In writing these equations we have used the following notations 

Mu = — , AMi = , Fi = ^ , 6.26 

rrie rUe rUe 

Pa ^ Pa ,„ „„N 

Pa = 7^ , Pa = 7^ , (6.27) 



with a = e, 7, v. 

As mentioned above, pe, Pe and p-y, p^, depend on the electron chemical potential. In 
the new numerical code we are writing, we use the following strategy. First of all, we 
have checked that the influence of (pe on the various quantities of interest is very weak, 
justifying the expansion at first order in (pe leading to Eq. (|6.18| ). By virtue of this fact. 



we have then substituted Eq. (|6.18D in the corrected expressions for pe, Pe, p-y, p-y (given 
by Eqs. (p. 53 ) - ( 2.54 ) plus Eqs. ( 5.18 ) - ( 5.21| )) and subsequently we have performed 
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a fit for these quantities as functions of z. To speed up the running of the BBN code 
we have finally used these fits in the evolution equations. The obtained values for the 
electron energy density and pressure are fitted (in the relevant temperature range) by the 
expressions 

Pe = 1.145 + 0.033981 z - 0.14543 + 0.025507 - (0.54168x 10~^) 

- (0.11263x10^^) - (0.29742x10"^) z^ + (0.38331x10"^) 

+ (0.45263x10"'^)/ + (0.19241x10"^) - (0.96597x 10"^°) 

- (0.19505x10"^°) - (0.14079x10"^^) , (6.28) 

Pe = (0.3786 + 0.019126 z - 0.063895 + 0.032085 z^ 

- 0.0048501 z^ - 0.00016611 z^ + 0.000082922 z^ + (7.9884x10"^) z'^ 

- (0.0619x10""^)^^ - (1.9568x10"^) z^ - (1.0921x 10"^) z^° 

+ (3.8564x10"^) z^i) e"°-™ . (6.29) 



The expression for z derivative of pe, entering in ( |6.23| ) and ( |6.24| ), are obtained from 

(H). 

Instead, in the considered temperature range, we have found that p^ takes values between 
0.6580 and 0.6573 while p^ varies from 0.2193 and 0.2187. For simplicity we have then 
taken the average values of 0.6577 and 0.2190 as the ones appropriate for p-y and p^ 
respectively. 

Finally, the most critical numerical part of the BBN code concerns the solution method 



of the set of differential equations ( |6.20D for Yi. In fact, at high temperatures, nuclear re- 
actions proceed in both forward and reverse directions almost equally rapidly, so that the 



right-hand-side of ( 6.20 ) is a small difference of large numbers, which causes severe numer- 
ical problems. To avoid this, we have used the Gear method of backward differentiation 
formulas for stiff problems, described for example in |93]. 



The implementation of the described code, using a Fortran Power Station, is currently 
under study, and represents the next step towards a precise determination of light element 
abundances in the Universe. 

The preliminary results obtained by our code for '^He mass fraction and D, '^Li abun- 



dances as functions of rj are shown in Figures 5.7, 6.8, 6.9 
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Figure 6.7: Our predictions for "^He mass fraction Yp versus rj. The solid line refers to 
Ni, = 3, while the dashed ones to N^, = 2 (below) and Ni, = 4 (above). 




Figure 6.8: Our predictions for D abundance Y2 versus ry (notations as in Figure 
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Figure 6.9: Our predictions for ^Li abundance I7 versus rj (notations as in Figure 
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Appendix A 

Non equilibrium distribution of 
species 



Decoupling of species from the background plasma is properly described by the evolution 
of the particles's phase space distribution function given by the Boltzmann 

equation. In this appendix we follow the approach of |Q] and restrict our analysis to an 
homogeneous and isotropic plasma, described by the Robertson- Walker metric. In this 
case / = f{\p\,t) or equivalently / = f{E,t) and the Boltzmann equation is 

where C[f] is the collisional operator for the processes creating and destroying the given 
particle specie. It is more convenient to use the number density n(t) rather than the 
distribution function, defined by 



(2^ 



= 9 I 7^fiE,t) , (A.2) 



g being the internal degrees of freedom. For this quantity the Boltzmann equation becomes 

with H = R/R. Let us focus on a given particle ip, interacting with the other particles 
in the heat bath through the reaction ^ + a + b + ... < — > i + j + for this process the 
collisional term can be written in the form 
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Pfe... - Pi - Pj...) X 



^ {\^\^ forward^ forward " {M]"^ f^^^j^^^^^ ^backward) , (A.4) 



with 



'^forward = fi,fafb-{l±fi){l±fj)- (A.5) 
^backward = fi fj ... (1 ± U) {1 ± fa) {I ± fb) ■■. (A.6) 

where +/- refers to bosons/fermions. The phase space density is given by 

dUr = dU^dUadUb...dUidUj , (A.7) 

while the squared matrix elements |Mp are averaged over initial and final spins and, 



assuming CP invariance, |M|2^„^^^^^ = |M|2^^^^^„^^ = |M|2. 

For illustrative purposes, we consider only situations in which Bose condensation or Fermi 
degeneracy can be neglected, so that we may safely disregard the stimulated emission or 



blocking factors in (A.4): (lit/a) ~ 1. In these approximations, the evolution equation 
( [A .31 ) simplifies to 

^ + 3Hn^= - J dUT{27^)U\p^ + pa + Pb...-p^-Pj...)JWx 

X iU fa fb... - fifj.:) . (A.9) 

Note that the 3H term accounts for the dilution effect of the expansion of the Universe, 
while the right hand side term accounts for the number changing ^/^ interactions. The 
dilution factor can be absorbed by considering the particle number in the comoving volume 
E? , proportional to the quantity 

y = ^ , (A.io) 

s 

where s oc is the entropy density. In fact, in this case the left hand side of ( [A.S| ) 
becomes simply 

-dT + '^"^ = 1^ ■ ^^-''^ 
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We are now in a position to describe the decoupling (i.e. the exit from chemical equilib- 
rium) of the given specie from the background plasma. Restricting ourselves to the case 
in which this specie is stable, the only number changing processes are of the type 

tp^ < — > XX . (A.12) 

With X we generically denote all the species into which ip can annihilate. These particles 
are, of course, in thermal equilibrium and hereafter we neglect, for simplicity, any chemical 
potential. Then we have 

=^ exp|-:?|^| (A.13) 
but, by means of energy conservation, the following condition 

/x - exp (- ^1 exp I- $U (A.14) 



holds, where with the index EQ we have indicated the equilibrium distributions. Hence 



the statistical factor in (A. 9) takes the form 



Uh - fxfx = Uh - f^'^f^'^ (A.15) 



and the Boltzmann equation for the present case becomes 



+ 3Hn^ = - < a^^^^j^v > - (nf^)^) (A.16) 



dt 

where < av > is the thermally averaged annihilation cross section times velocity: 



Iq- J dUr (27r)^ 6^{p^ + - Px - P^) 



< f^v^^xx^ > - f EQ 

X Pwf exp|-:|^}exp|-:^| . (A.17) 

We stress that the rate of change of ip is proportional to the annihilation rate, while 
it tends to zero when the equilibrium is approaching since creation processes balance 
destruction ones. 

Another impressive form of ( [A.16 ) is obtained by considering the Boltzmann equation for 



the quantity Y and rewriting the t dependence in terms of the temperature or, better, in 
terms of the adimensional parameter z = m^/T: 

2 



z dY 



Yeq dz 

104 




1 (A.18) 



r = UEQ <av> . (A.19) 

Prom these we immediately see that the rate of change of tp per comoving volume is set 
by the factor F/H: when this is less than order unity, the relative change in the number 

of ip in the comoving volume becomes small, annihilations freeze out and the abundance 
of il) "freezes in" . 
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Appendix B 

Neutrino-photon temperature 
relation 



When neutrinos decouple (at TJ^ around 1 MeV) from the background primordial plasma, 
their thermal distribution retains the equilibrium Fermi-Dirac form with a temperature 
Ty scaling (approximately) as the inverse of the cosmic scale factor R. As discussed in 



1.4.1, their entropy Sy in the comoving volume is conserved also after decoupling 

5,(ri^) = S,{Ty) . (B.l) 

On the other side, photons still remain in equilibrium with the pairs until these disap- 
pear when become non relativistic (around T ~ rrie). However, since the total entropy in 
the comoving volume is conserved as well as neutrino entropy, it follows that the entropy 
Si of the species still in equilibrium (/ = 6^,7) is conserved too 

Si{Td) = Si{T) . (B.2) 

Note that, until neutrinos decouple, = T; consequently T^P = T^. According to Eqs. 



( |L77|) and (IlIqI) , (|L50|) , the entropy of a specie is explicitly given by 



5, = ^I{x,)R^Tf , (B.3) 

/(X.) = r,^(^^ + 2x.,)V^(V + 8x,., + 3xf) 
Jo e^t+y ± 1 

where Xi = mi/Ti, y = {Ei — mi)/Ti. Here, the effective (T dependent) degrees of freedom 
are parametrized by gi(T) = {gi/GTr"^) I{xi). Substituting ( B.3 ) in ( |B.1D , ( |B.2| ) and dividing 



side by side these two equations, the following relation between and T is obtained [92| 

'Ty\' _ 9u{Td) 9i{T) 



Tj gi{TD) gAT) 
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Figure B.l: Neutrino - photon temperature difference, 
or, more explicitly, 



% f 2Tr^ + 151{xe 



27r4 + 151{x^] 



(B.6) 



In the limit (considered in |l.4.l| ) in which gi{T) ~ 2, giiTo) ^ 11/2 {guiTo) = gu{T)), 



from (^) the relation ( |1.119|) follows. Instead (^ ) gives the generalization for Ty as 



function of T we looked for. In Figure we plot the relative difference in temperature 
between neutrinos and photons; of course, before decoupling the two temperature coincide 
(the plotted curve is evaluated from ( |B.6| ) adopting the value 2.3 MeV for the decoupling 
temperature; results practically unchange by shifting this temperature in the range 2 
3 MeV). 
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Appendix C 

Nuclear reaction network 



In the following Tables C.l 



we show the complete nuclear reaction network relevant 
for the calculations of the primordial abundances of light elements 62 , ^ . 
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Table C.l: Nuclear /3-decay reactions. 









+ Fe 


+ -^He 






e~ 


+ Fe 


+ 24i^e 






e~ 


+ !7e 


+ 






e~ 


+ !7g 


+ 14 AT 






e+ 


+ l^e 


+ 24ife 






e+ 


+ t'e 


+ lis 






e+ 




+ 


13iV 




e+ 


+ Ve 


+ 13c 


14(5 




e+ 


+ l^e 


+ i^iV 


15^ 




e+ 


+ I^e 


+ i^iv 



Table C.2: Nuclear neutron-photon reactions. 
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(n 


7) 


D 


D 


(n 


7) 




^He 


(n 


7) 


^He 


^Li 


(n 


7) 


''Li 


''Li 


(n 


7) 


""Li 


10^ 


(n 


7) 


115 


lis 


(n 


7) 


125 


12c 


(n 


7) 


13c 


i^C 


(n 


7) 


14 c 


14iV 


(n, 


7) 


15 AT 



Table C.3: Nuclear neutron-proton and proton-neutron reactions. 
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lie 
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137V 


{n,p) 
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Mat 


{■n,p) 
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12^ 


{p, n) 


12c 



109 



Table C.4: Nuclear neutron-a and a-neutron reactions. 
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Table C.5: Nuclear proton-photon reactions. 
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Table C.6: Nuclear proton-a and a-proton reactions. 





(P, 


a) 




'Li 


{P, 


a) 


^He 


^Be 


{P, 


a) 


^Li 


10^ 


{P, 


a) 


'Be 




{P, 


a) 






{P, 


a) 


12^. 




(a 


P) 




10^ 


(a 


P) 


13^ 




(a 


P) 






(a 


P) 






(a 


P) 




13 AT 


{a,p) 


1%Q 



Table C.7: Nuclear a-photon reactions. 
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Table C.8: Nuclear deuterium-nucleon reactions. 
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Table C.9: Nuclear three-body reactions. 
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Appendix D 

Nuclide mass excess 



In this appendix we report the measured mass excess Amj for the 26 nuchdes used in 



BBN calculations |35]. 
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Table D.l: Nuclear mass excess Amj in KeV (referred to ^^C). 



n 8071.323 ± 0.002 


8667.984 ± 0.420 


p 7288.969 ± 0.001 


^^C 10650.531 ± 0.952 


D 13135.720 ± 0.001 


^'^B 13368.901 ± 1.400 


■^H 14949.794 ± 0.001 


0.0 


'■^He 14931.204 ± 0.001 


^^Af 17338.083 ± 1.000 


^i^e 2424.911 ± 0.001 


^^C 3125.011 ± 0.001 


^Li 14086.312 ± 0.475 


^•^N 5345.456 ± 0.270 


■^Li 14907.673 ± 0.473 


^'^C 3019.892 ± 0.004 


''^e 15769.489 ± 0.472 


14 AT 2863.417 ± 0.001 


«Li 20946.195 ± 0.488 


i^O 8006.456 ± 0.075 


«S 22921.002 ± 1.107 


^^N 101.438 ± 0.001 


^Be 40818.362 ± 62.471 


2855.388 ± 0.503 


12050.761 ± 0.370 


i^^O -4736.998 ± 0.001 



114 



Appendix E 



Finite nucleon mass corrected 
matrix element 

In this appendix we report the result for the spin summed squared modulus of Eq. ( [4.lD 
J2\M\^ = 16 ( Ca^ {Ml M| + 2 Ml Ml + M| M| + 2 Ma M4 + 2 M2 M| M4 + 

spins 

+ Ml Ml + 2 M| M| + M| M| - 2MIs - 2 M| s + 

- 2MIs - 2MIs + 2s^ - Mlt - M^t + 

- M| t - 2 M2 M4 i - M| t + 2 s t + t^) + 

+ (Mf M| + 2 M| + M| M| - 2 M2 M4 - 2 M2 M| M4 + 

+ m| + 2 M| M| + M| M| - 2 s - 2 M| s + 

- 2 M| s - 2 M| s + 2 - f - M| t + 

- M| t + 2 M2 Mit - Mjt + 2st + t^) + 

+ 2CAfps{- Ml M2 Ml + M| M| + M2 M| + M4 - Ml M| M4 + 

- M| M| M4 - Ml M2 Ml + Ml + M2S - M2 Ml s + 

- M4S + Ml M4S - M2Mlt - Ml M4 1) + 

- 2Ca — {Ml M| - M| M| + M| M4 - M| M| M4 - Ml M2 Ml + 

Mn 

+ M2 Ml Ml - Ml Ml + Ml Ml - Ml M2t - M^t + 

- M2 Ml t - Ml Mit - Ml Mit - MlM^t - M2Mlt + 

- Mlt + 2M2st + 2MiSt + M2t^ + Mif) + 

+ 2Cvh{- Ml M2 Ml + Ml Ml + M2 M| - Mf M4 + Ml Ml M4 + 
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+ M| M| M4 - Ml M2 Ml - Ml Ml + Ml M2S - M2 M| s + 

+ Ml M4S - M^ M4S - M2M^t + Ml M4 1) + 

+ 2CvCa{- Ml M| + M| M| + M| - M| M| + t + 

+ M'^t + M^t + Mlt - 2st - + 

+ 2 — (Mf M| + Ml M2 Ml - M2 M| - Mf M4 - Ml M| M4 + 
Mm 

+ m| M4 - M2 Ml Ml + M| M| - M2S + M2 Ml s + 

+ Ml M4S - Ml M4S - Ml M2t - M^t + Ml M4 f + 

- M| M4 1 + M2 M| t - Mft + M2t^ + M4 1^)) (E.l) 

In the above formula, neutrino mass has not been set to zero to give a general expression 
holding for all reactions. For the two-body scattering processes in (p.2|), ( p. 3D it is 

s = {pi + p2f = (p3 + pif , 
t = {Pi - Psf = {P2 - Pif ■ 

As in section (p.l|), pi and p^ are the initial and final lepton momentum while p2 and p^ 
are the initial and final nucleon momentum, respectively. Instead for n ^ pe" V^, it is 

s = {-Pi+P2f = {p^ + Pif , 
t = {-Pi-P^f = {-P2+PAf ■ 

while lov pe~ Vf. — > n it is 

S = {Pl + P2f' = (-P3 + Pif 
t = (pi + p^f = {-P2 + Pif 
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Appendix F 



Finite temperature and density 
QED 



F.l The Real Time Formalism 

At finite temperatures, physical processes take place in a heat bath consisting of a back- 
ground of a plasma of particles and antiparticles. All interaction processes must be there- 
fore considered as occurring in a thermal "vacuum" state, which does not coincide with 
the usual Poincare invariant vacuum of zero temperature quantum field theory. Manifest 
Lorentz invariance therefore breaks down due to the choice of the preferred frame of the 
bath, and the zero temperature renormalization prescription cannot be applied straight- 
forwardly. The covariant formalism that enables to consider the interactions of a given 
particle with the surrounding medium is the Finite Temperature and Density Quantum 
Field Theory | p6| , |8^. Here we adopt the Real Time formulation of this theory, in 
which the Feynman rules for the vertices are identical to the corresponding ones in the 
vacuum. The effect of the temperature and of the density is taken into account in the 
expressions of the free particle propagators. For fermions and bosons we have respectively 






where 



r(p) = 2Ti 5{p^ 



rm?) {9{p-u)n{p) + 9{—p-u)n{p)) 



(F.3) 
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and 

(P = 1/T) are the Fermi-Dirac and Bose-Einstein distribution functions (np and ub are 
the distribution for the antiparticles, obtained replacing fi with — /u). Note that in a 
plasma, another 4-vector must be considered, the 4-velocity of the medium, which in 
its rest frame is given by u'^ = (1,0). 

The additional contributions in the above formulas (denoted with Tp^s) pick up real par- 
ticles through the mass shell (5-functions and are proportional to the particle densities in 
the thermal bath. They take into account the role played by the medium in the single 
particle propagation. 

The heat bath we are interested in is the primordial plasma at the epoch of BBN, i.e. at 
temperature around 1 MeV. At these temperatures the and Z'^ gauge bosons degrees 
of freedom are not excited, so we can safely neglect the temperature dependent term Tb{p) 
in D^u for these bosons. 

Our primordial plasma consists of nucleon, pairs, photons and neutrinos. The temper- 
ature dependent term in the propagator of nucleons would be suppressed by a Boltzmann 
factor smaller than exp(— 100), compared with the ones for the other particles, in the 
temperature range relevant for BBN, so it will be neglected in the following as well. 
The most important finite temperature effects induced by the remaining particles in the 
heat bath come from the QED interactions of and photons, which are of order a. Neu- 
trinos only interact weakly, thus their contribution is at least of order Gp. 
In this appendix we calculate the self-energy and wavefunction renormalization at 
finite temperature and photon self-energy, which are relevant for the thermal radiative 
corrections considered in chapter H. 

F.2 Electron self-energy 



The Feynman diagram for electron self-energy S is reported in Figure F.l; it gives 



= I ^{^er)iSp{p-k){ier)iD,,{k) (F.6) 
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Figure F.l: Self-energy diagram for an electron in a plasma. 

where Splp—k) and D^y{k) are the (RTF) Fermi and Bose propagators in a plasma, which 
are reported in ( |F.1| ), ( [F.2| ). Here we are only interested in the temperature dependent part 
St of the self-energy; moreover, we neglect possible damping effects (which are suppressed 
with respect to coherent ones) and then consider only its real part: 

^(^t) = + = / (0 ((i>- ^ - 2m.) _\fl + 

TF{k) 
{P - kf 



- 2m,) j::^^] (F.7) 



From Lorentz invariance, it follows that the self-energy has the general form 

JJSt = ai) + h^ + c (F.8) 

where a,b,c are functions of the only two invariants ,p-u (obviously = 1). In the 

following, when necessary, we shall adopt the plasma rest frame, where = (1,0). We 
then have 

a = -\{Tp- poTu) (F.9) 
P 

b = ((Po -P^Tu- poTp) (F.IO) 

c = ^rrJRSr (F.ll) 

where 

Tp = ^Tri^^J^T) (F.12) 
Tu = ^TriA^^T) (F.13) 
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For the Bose parts, after some calculations we obtain (p = |p|, /c = |k|) 



p 



J dkkB{k) (^-4 + {h{p,,P.k) + h{-po,p,k))^ (F.14) 

= -^J dkB{k) {{po - k)h{po,p,k) +P0 + k)h{-po,p,k)) (F.15) 

CB = / dkB{k) {Ii{po,P,k) + hi-po,p,k)) (F.16) 

vr J 

where tj = Pq — p'^ — ml. Instead for the Fermi parts we get: 

= -^Jdk^F{E) {hipo,P,k) + l2{-po,P,k)) (F.17) 

7f = -^J dkk''F{E) {h{po,p,k) - h{-po,p,k)) (F.18) 

CP = ^ fdk'^FiE) {l3ipo,p,k) + h{-po,p,k)) (F.19) 

TT J E 

with E = \/p2 + rrig. The quantities /j are defined by ^ 

/i(P0,P,A:) = ^ = - In " t ^ ^1 (F.20) 



2{pqE — pkx) 



rj + 2mg — 2poE + 2pA;x 

,2 1^ I 



h{po,P,k) 



r] + 2m^ |r? + 2mi - 2poE + 2pfc| ^ 
2pA: |r? + 2m2 - 2poE - 2pk\ ^ ' ' 

1 2dx 

-1 ?7 + 2mg — 2po-E' + 2pkx 

pk \r] + 2m2 - 2pqE - 2pk\ 
In the above formulas we have treated pq and p as two independent variables. Since the 
expressions in ( F.14| )-( |F.19| ) are already of order a, in these we can substitute the vacuum 



dispersion relation ~ + "^e? thus obtaining the following results for the coefficients 
a, 6, c 

' -1 + — In + 



V 2p u) — p 



^We have chosen a reference frame in which the electron momentum hes along the z axis [x — cos 
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a 



— ^ a; 1 - In + 



a me 
irp 



J dk^F{E) {In A - InB) , (F.25) 

^ ^ + ml +pk ^ ^ ujE - ml + pk 

LoE + ml — pk ' ujE — ml — pk ' 

with UJ = \/p^ + ml. 

From the calculated electron self-energy, we can now proceed to evaluate the mass shift 

correction. 

If we denote with ip the wavefunction of the electron in the plasma, it satisfies the 
correct Dirac equation 

(i» - me - S) ^ = , (F.27) 
with S in ( [F.8| ). We can rewrite this in the form 

{p - m) = , (F.28) 

with 

p^ = (1 - a)pf, - bUf, (F.29) 

m = me + c (F.30) 



By left-multiplying Eq. (|F.28D hy J) + fh we then deduce the correct dispersion relation 



for an electron in a medium 

p2 _ ^2 ^ , (F.31) 

or 

Po - P' - 1 Po + \^ ' = . F.32 

1 — a (1 — ay 

At first order in a we finally obtain the energy shift 

Po - \ p^ + ml + [a^ + h + c— \ = UJ + ^l , (F.33) 

* \ UJ UJ 



where a, 6, c are a, 6, c evaluated in pQ = uj. 
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F.3 Wavefunction renormalization 



In general, the Dirac spinors corresponding to a particle propagating in a plasma are 
different from the ones travelling in vacuum. This results in effective temperature de- 
pendent energy projection operators which replace the usual vacuum ones in the 
computation of spin summed squared amplitudes 

= - , (F.34) 

and similarly for the negative energy projector. Thus, we first have to calculate the 
effective energy projection operators and then evaluate the corresponding corrections for 
the decay rates for the reaction in (^), (p.2p, ( p.3P . 

An electron propagating through a plasma is described by the equation of motion in 
( [F.27 ), so that the field propagator is given by 



p — rrif. — Tj P — m — fh? 

We identify the particle states as corresponding to the energy poles in the propagator; the 
wavefunction renormalization can then be read off by evaluating the residue of G at the 
pole. For definiteness, let us focus on the positive energy projection operator; a similar 
procedure will hold for the negative energy one. If we expand G in ( |F.35| ) around the 
positive energy pole, say ojr, we obtain 

-1 / . o . \ -1 



— m^j = — rrig — 2ap^ — 2bp-u — 2crnj 

~ [pI- uj^ - 2{a + a' {po - ujr)) (pi - p^) + 

+ 2 (b + y {po - ujRjj p-u - 2m {c + c' (po - ^i?,))} ^ 

~ |po - uj"^ - 2apl - 2 a' uj^ {po - ujr) + 

+ 2p^ (a + a {po - ujr)) - 2b p-u - 2y p-u {po - ujr) + 

- 2m {c + c {po - ujr))}~^ . 

Here we have used the "hat" to denote that a given quantity is evaluated at po = ujr, while 
a "prime" is used to indicate differentiation with respect to pq. Note that throughout the 
entire calculation we retain only terms which are of the first order in a. Substituting the 
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relation bp-u ~ bp-u + buQ (pQ — loji) in the above formula, we get 



(^p^ - m^^ ^ = |(l-2a) pI - (^LJ^ + 2aml + 2bp-u + 2cm 
- {po -iOR) 2b uo + 2a' ml + 2b' p-u + 2c'm | 



+ 



By requiring that G (and then (p^ — m^) ^) has a pole for po = ujr we deduce that 



uj\ = uJ^ + 2am1 + 2bp-u + 2cm , 



(F.36) 



which exactly corresponds to the energy shift in ( [F.33 ). Around po = ujh we then obtain 



/ ~2 ~ 2 

[p — m 



-1 



1 



— \ l + 2a + - UQ + a' ^ + b' — + d — 

Pq — UJji ZLOr \ UJ UJ UJ UJ 



(F.37) 

thus the expansion around the positive energy pole pQ = ujr of the propagator G in ( |F.35| ) 
finally gives 



G 



1 



yp-u ^,me\ i> + mR 



Po - 
b 



+ 



1 + aH uq + a! — + b he — 

U) U U UJ 



UJ 



2ujR 



2uj \ me 

Here, we also report an alternative writing of the same expression 



(F.38) 



G 



1 



Po - 
b 

2J^ 



l + a+-^(no-l)+a'^ + 6'^ + c'^ 

UJ UJ UJ UJ 



2UJR 



+ 



(a; /i ■ 



+ 



(F.39) 



The quantity mR is the effective mass corresponding to the effective energy ujr in ( |F.33 ) 



mR ~ me + dm ~ m^ + am^ + b + c 

me 



(F.40) 



In the particularly relevant case in which = (1, 0) (plasma rest frame) we obtain 



G 



1 



mt 



, me \ p + mR 



Po - UJR 



1 + a + a' ^ + b' + c' — , 

UJ UJ I 2ujR 



^ I P-7 + — 

2uo^ \ mp 



(F.41) 

The positive energy projection operator is then simply the residue of G at the correspond- 
ing pole; for example, in the plasma rest frame we have 



A+ = (1 + A) 



i> + mR 
2ujR 



2u;2 



P-7 + — 

me 



(F.42) 
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with 



A = a + a' ^ + 6' + c' ^ . (F.43) 

UJ LO 



For the negative energy projector, one has simply to make the substitutions w —lo, 
ujR — > —ujR-, p — > — p in all the quantities appearing in ( F.42| ), ( F.43 ) (in particular these 
substitutions have to be made also in the expressions for a, a', ...). 

Some considerations are now in order. We mainly note that differently from the usual 
vacuum case, wavefunction renormalization at finite temperature is not simply given by 
a multiplicative renormalization factor (in the vacuum case A/j = (1 + Z2) Aq), but also 
introduces an additional term proportional to h and with a different spinorial structure. 
This is a pure matter effect, which is given by the presence of the medium 4- velocity 
Ufj_. This dependence is, however, also present in the multiplicative renormalization factor 
(the first term between {...} in ( F.38| ) or ( F.39| ) or in A in ( F.42| )); we point out that, 
in the plasma rest frame, the peculiar medium dependence of A only occurs through the 
energy derivative of 5, differently to what happens in the general case in which also a 
term proportional to h is present (see the first term in ( [F.39D ). One can easily check that 
in the limiting case of no medium, we recover the usual expression for the renormalized 
projection operator. We note, however, that this is not true |Q] for the result quoted 
in [^2|, although an approach similar to ours is used. Other different approaches to the 
problem of wavefunction renormalization have been proposed in literature ]7^-[pl|], but 
the final results striking differ. For example, in WM the authors start introducing finite 



temperature spinors chosen as the solution of the (nonlinear) Dirac equation ( F.27 ) whose 
corresponding creation and annihilation operators are assumed to satisfy ordinary, zero 
temperature, anticommutation relations. Expanding the propagator in terms of these 
spinors they obtain a wavefunction renormalization factor which is only multiplicative. In 
particular, this multiplicative factor is exactly the same as our factor (1 + A) in (p\42D , but 
the additional term proportional to h is absent. The essential difference with our approach 



is the assumption made in [74| on the canonical spinor basis to be used; here, however, we 
do not make any hypothesis on the renormalized field to be used, but simply recover the 
particle content and the corresponding projector operators from the poles of G and their 
residues, respectively. Comparative analysis of different approaches have been carried out 
inRefs. Isl. 



124 



p 

Figure F.2: Self-energy diagram for a photon in a plasma. 



F.4 Photon self-energy 

The Feynman diagrams for photon self-energy II^i, is shown in Figure [F.2| ; it gives 

iU^, = - J -0-^Tr{{ier)iSF{p) {ieY)iSF{p-m , (F.44) 

where again Sf{p) is the electron propagator in the plasma given by ( [F.l| ). As for the 
electron self-energy, we are interested only in the real part of H^i, which, after some algebra, 
is given by 

2 f d^p 1 (2p^py - Y>iiK - Pukf, + p-kgi_,y 



, '^P^lPu + VtiK + Puk^, - p-kg^u ^ 

+ W^^k )K + n,) (F.45) 



where E = \/'p^ + m'j. It is now convenient to express Hf^p in terms of form factors by 
using Lorentz and gauge invar iance, which give the following constraints 

fe'^n^, = rn^, = . (F.46) 

In presence of a medium, II^i, may be expressed in terms of the tensor g^i^ and the 4- vectors 
k^ and u^. Then, noting that 

g^lu = g^lu - , (F.47) 



and 



(F.48) 
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with Ufj^ = gi_iuu'^ , are the only two tensors orthogonal to A;^ that can be built with g^,j, 
and u^, the most general form of 11^;^ (assuming parity conservation) results 

= IlTR^,u + HlQ^u , (F.49) 

with 

Ruu — Eluu T'n ) Qau — T'n • (F.50) 

The quantities IIt, ITj^ are scalar functions of the only two invariants k'^, k-u and in terms 
of li^i, are given by 

= g'^^n^. , (F.51) 

2Ilt = R^^'Ii^. ■ (F.52) 



It is now simple to explicitly evaluate these form factors from ( F.45| ); at first order in a, 



substituting the vacuum dispersion relation — = in the expression ( [F.45| ), after 
some algebra we obtain 

^{Hl) ^ O(a^) , (F.53) 



2a p 



2 



^{Ut) — dp ^ / ,^ (ne + ne) + 0{a') . (F.54) 
From the equation of motion for the 4-vector potential field 

{-k'g^,, + n^.) A" , (F.55) 
which more usefully can be written more usefully in the form 

(A;2 - Ut) R,,u + {k^ - Ul) Q,,u = , (F.56) 
we can deduce the perturbed dispersion relation of the transverse modes, which is 

uj^ -k^ = Ut , (F.57) 
with given, at first order in q, in ( [F.54| ). 
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